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CHAPTER  1 
INTRODUCTION 


It  is  the  purpose  of  this  monograph  to  present  a  critical  review  of  the  ability 
of  the  engineering  profession  to  anticipate  the  dynamic  environment  of  a  flight 
vehicle,  to  predict  the  vehicle  response  to  this  environment,  and  to  simulate  the 
expected  environment  in  the  laboratory.  The  ability  to  do  these  things  has 
progressed  significantly  in  the  last  few  years,  although  engineering  practice  has 
not  completely  kept  up  with  the  rapidly  advancing  technology.  It  is  hoped  that 
this  monograph  will  aid  in  spreading  this  new  information  to  a  larger  group  of 
environmental  engineers. 

Many  of  the  methods  described  in  the  following  chapters  will  rely  on  argu¬ 
ments  motivated  on  physical  grounds  as  well  as  on  mathematical  models  of  the 
environment  and  the  structure.  Environmental  engineers  have  generally  taken  the 
attitude  that  predictions  of  response  must  be  made,  even  in  the  face  of  imperfect 
data  and  methods.  For  this  reason  some  empiricism  is  evident  in  the  arguments. 
Empiricism  is  frequently  useful  because  of  the  economics  of  effort  that  can  re¬ 
sult  from  simple  rule-of-thumb  procedures. 

Generally,  the  studies  that  we  report  on  are  applicable  to  sustained  loads  that 
are  random  in  time  but  may  or  may  not  be  spatially  homogeneous.  We  exclude 
by  implication,  therefore,  such  loads  as  pyrotechnic  shock  and  combustion  or 
structural-propulsion  interaction  instabilities.  In  certain  cases,  it  is  possible  to 
treat  these  loads  by  statistical  methods,  but  the  appropriate  conditions  require 
experience  and  judgment  for  their  application. 

The  following  paragraphs  describe  briefly  the  subject  matter  of  the  book.  In 
that  sense  they  serve  as  an  outline.  Their  primary  purpose,  however,  is  to  make 
explicit  the  interconnecting  ideas  and  the  logical  relationships  within  the  ma¬ 
terial  to  be  covered. 

1.1  Prediction  of  the  Environment 

In  the  totality  of  the  environmental  problem  comprising  prediction  of  loads, 
calculation  of  response,  and  design  of  simulation,  the  loads  prediction  problem 
is  by  far  the  least  advanced.  In  part,  this  is  due  to  our  inability  to  predict  How 
behavior  from  the  Navicr-Stokcs  equation  [1].*  It  is  also  due  to  the  multitude 
of  different  but  significant  loads-acoustic,  turbulent  boundary  layer,  wake 


*  References  are  cited  in  the  Bibliography,  appendix  B. 
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impingement,  oscillating  shocks,  separated  flows,  etc.  [2]  -that  are  experienced 
by  the  vehicle.  These  loads  appear  on  a  vehicle  at  different  places  during  different 
periods  of  its  trajectory  and  their  intensity  and  excitation  characteristics  depend 
on  both  trajectory  and  geometric  parameters  [2,  Sect.  2] . 

There  are  three  ways  that  information  on  loads  is  obtained:  through  theo¬ 
retical  computation,  laboratory  studies  on  models,  and  flight  measurements  on 
vehicles.  As  a  practical  matter,  a  particular  estimate  is  often  a  combination  of  all 
three  methods.  Theoretical  and  laboratory  studies  suggest  general  features  such 
as  overall  levels,  spectral  shapes  and  scaling  laws,  while  flight  measurements  on 
similar  vehicles  serve  to  corroborate  these  studies  and  to  indicate  the  kinds  of 
loads  that  may  be  present  [3] . 

The  equations  governing  structural  motion  require  considerable  detail  in 
specifying  the  applied  loads,  either  in  terms  of  pressure  as  a  function  of  time  at 
all  positions  on  the  surface,  or  by  statistics  of  the  pressure  such  as  the  space-time 
covariance  [4] .  Since  only  a  limited  number  of  measurements  can  be  practicably 
obtained  in  laboratory  or  flight  studies,  it  is  necessary  to  develop  conceptual 
models  of  the  environment.  These  are  fitted  to  the  expected  environment  by 
laboratory  and  field  measurements,  but  they  always  imply  more  (and  less)  detail 
than  can  be  justified  from  the  data  [5].  We  shall  explore  this  point  more  fully 
in  the  discussions  to  come. 

In  summary,  it  is  possible  at  present  to  give  a  description  of  the  dynamic  en¬ 
vironment  for  a  particular  vehicle  in  a  given  trajectory.  This  description  will  have 
its  weak  points,  and  much  work  is  required  to  improve  it,  but  it  can  serve  as  a 
basis  for  response  calculations.  We  turn  to  that  subject  now. 

1.2  Calculation  of  Response 

An  aerospace  flight  vehicle  is  a  very  complicated  structure,  consisting  of  the 
outer  dixiii,  ii4iiivo  and  stingc-is,  and  internal  structure  and  acoustic  spaces.  In 
quite  a  real  sense,  it  defies  description,  except  perhaps  through  the  library  of 
technical  drawings  used  in  the  specification  of  its  construction.  Even  if  one  had 
an  exact  mechanical  description  of  the  structure,  however,  with  every  rivet  and 
weld  exactly  positioned  and  defined,  and  a  magnificent  computer  able  to  comain 
all  of  the  structural  and  environmental  data,  then  one  could  still  be  dissatisfied 
with  solutions  of  response  obtained. 

The  development  of  scientific  or  engineering  understanding  of  a  system  is 
closely  related  to  one’s  ability  to  reduce  tire  system  to  simpler  elements.  To  the 
extent  that  one  can  identify  certain  response  characteristics  as  due  to  cylinder 
behavior,  or  panel  modes,  or  fluid  loading,  it  is  possible  to  predict  how  the 
response  might  be  changed  by  changing  appropriate  properties  of  the  structure 
(2,  Sect.  3] . 

In  the  following  discussions,  the  methods  for  computing  response  are  clas¬ 
sified  in  accordance  with  the  way  that  the  stiucture  is  modeled.  A  given  structure 
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can  be  modeled  in  several  ways,  depending  on  the  frequency  of  response  and  the 
kind  of  excitation.  It  is  the  engineer’s  job  to  pick  the  proper  structural  and  en¬ 
vironmental  model  for  the  calculation  of  response. 

By  and  large,  the  models  ignore  effects  of  nonlinearity  in  response  [6] .  A 
major  distinction  can  also  be  drawn  between  models  that  are  infinite  in  extent, 
and  those  that  are  finite.  The  former  generally  require  calculation  of  traveling 
wave  response,  while  the  response  of  the  latter  is  expressed  as  vibration  of  nor¬ 
mal  modes.  In  many  cases  it  is  possible  to  obtain  equivalent  results  from  these 
apparently  quite  different  systems. 

The  major  purpose  of  the  structural  model  is  to  predict  the  excitation  of  that 
part  of  the  structure  directly  exposed  to  the  environment.  Frequently,  however,  < 

one  is  also  concerned  about  mechanical  and  acoustic  energy  transmitted  to  in¬ 
terior  portions  of  the  vehicle  [7j.  Structural  models  for  transmission  generally 
require  a  more  subtle  design  than  those  designed  for  response  only.  The  dif¬ 
ference  lies  in  the  role  of  nonresonant  motion,  which  must  be  adequately  rep¬ 
resented  in  tiansmission  studies  [8].  We  shall  discuss  a  model  for  both  sound 
and  vibration  transmission  in  chapter  3. 

Finally,  we  must  recognize  the  role  of  laboratory  and  field  studies  in  response 
estimation.  In  fact,  much  of  the  cuuent  estimation  of  environmental  response 
is  based  on  purely  empirical  understanding.  Only  recently  have  analytical  methods 
become  capable  of  providing  a  reasonable  alternative  (or  complement)  to  the 
purely  empirical  procedures  [2]. 

1.3  Simulation  of  Environmental  Loads 

The  simulation  of  a  load  is  the  replacement  of  that  load  by  another  designed 
to  be  equivalent  to  it  in  some  sense.  The  simulation  may  attempt  to  copy  the 
service  load;  e.g.,  a  turbulent  boundary  layer  (TBL)  on  a  flight  vehicle  may  be 
simulated  by  a  TBL  in  a  wind  tunnel  model  [9].  The  simulation  may  replace 
one  load  with  another:  acoustic-induced  vibration  in  a  service  vehicle  may  be 
simulated  by  shaker-induced  vibration  in  the  laboratory  (10].  Whatever  the 
method,  it  is  desirable  that  a  good  understanding  be  achieved  of  the  phenomenon 
that  is  being  simulated. 

Suppose  we  arc  interested  in  the  noise  generated  within  a  spacecraft  shroud 
during  the  high-0  portion  of  its  trajectory.  Let  us  further  suppose  that  the  shroud 
vibration  had  been  monitored  during  this  portion  of  the  flight.  If  the  same 
vibration  levels  are  induced  in  the  shroud  by  placing  it  in  an  acoustic  chamber 
(and  thereby  assuming  an  equivalence  between  aerodynamic  and  acoustic  re¬ 
sponse),  then  it  may  happen  that  the  sound  levels  within  the  shroud  are  not 
correctly  reproduced  (11]. 

When  service  acoustic  loads  arc  simulated  by  sound  fields  in  the  laboratory, 
then  relatively  simple  scaling  laws  are  operative.  The  scaling  laws  for  some 
aerodynamic  loads  such  as  separated  Hows  and  oscillating  shocks  are  more  dif¬ 
ficult  to  demonstrate.  In  these  cases,  it  is  necessary  to  rely  heavily  on  data 
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taken  over  a  range  of  scale  factors  in  order  to  demonstrate  that  extrapolation  is 
reasonable.  Examples  of  such  extrapolation  will  be  given  [12] . 

The  principal  point  to  be  made  is  that  simulation  must  be  carried  out  with 
great  care  and  with  understanding  of  its  purpose  and  limitations.  Both  environ¬ 
mental  and  structural  parameters  will  affect  the  success  of  a  simulation.  Scaled 
attempts  to  replace  a  flight  environment  with  a  wind  tunnel  experiment  may  be 
quite  misleading,  while  another  test  of  panel  response  to  an  oscillating  shock 
might  be  quite  well  simulated  with  a  mechanical  shaker.  It  is  only  required  that 
one  understand  the  system  sufficiently  well.  How  well  the  systems  are  understood 
that  we  are  considering  will  be  revealed  in  the  following. 


CHAPTER  2 

FLIGHT  VEHICLE  NOISE  ENVIRONMENT 


In  this  chapter,  we  consider  several  aspects  of  the  noise  environment  of  flight 
vehicles:  what  the  major  types  of  environment  are,  how  they  are  measured,  when 
they  may  be  expected,  and  how  we  are  to  estimate  the  intensity  of  each  kind  of 
environment. 

The  kinds  of  environment  are  first  separated  into  acoustic  and  aerodynamic. 
By  acoustic,  we  mean  pressure  fluctuations  that  propagate  with  the  speed  of 
sound  and  are  related  to  velocity  fluctuations  by  the  formulas  of  ordinary 
acoustic  theory.  Aerodynamic  noise  is  a  broader  term,  encompassing  pressure 
loads  due  to  turbulence  and  shock  waves.  In  the  present  context,  it  includes 
turbulent  boundary  layers,  separated  flow,  wake  flow,  base  pressure  fluctuations, 
and  oscillating  shocks. 

The  reason  for  an  environmental  estimate  is  primarily  response  estimation  or 
simulation.  Calculations  of  response  require  considerable  information  regarding 
the  temporal  and  spatial  characteristics  of  the  loads-more  information  than  is 
generally  available  from  flight  measurements.  The  needed  information  is  usually 
inferred  by  fitting  the  data  to  an  assumed  model  of  the  load,  a  model  that  may 
be  established  from  laboratory  studies.  Insofar  as  possible,  we  shall  try  to  indi¬ 
cate  in  the  following  just  how  well  the  various  environmental  models  are  estab¬ 
lished  by  laboratory  and  field  data. 

2.1  Acoustic  Environment 

The  acoustic  environment  of  a  flight  vehicle  is  generated  by  the  engines,  by 
flow  over  the  vehicle  structure,  and  by  moving  parts  of  the  vehicle,  such  as 
propellers.  The  sound  generated  by  propellers,  rotating  helicopter  blades,  and 
the  compressor  blades  associated  with  bypass  jet  engines  is  now  reasonably  well 
understood.  Noise  from  this  source  is  generally  of  concern  from  the  community 
reaction  viewpoint  and  is  not  important  as  a  structural  load  [13].  The  noise 
caused  by  flow  over  the  structure  is  of  great  interest  in  terms  of  cabin  noise 
levels  (which  we  cover  in  sections  on  response  to  aerodynamic  noise).  In  at  least 
one  rather  special  case,  the  anticipated  sound  levels  radiated  onto  a  laminar  wing 
by  turbulence  is  of  interest  since  the  laminar  boundary  layer  on  such  an  airfoil 
might  be  quite  sensitive  to  noise  [14].  For  this  report,  we  shall  concentrate 
on  the  broadband  sound  radiated  by  subsonic  and  supersonic  jet  and  rocket 
engines. 
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Subsonic  Jet  Noise 

The  primary  interest  in  the  jet  noise  of  aircraft  has  been  the  radiation  of 
sound  away  from  the  vehicle  [15]  because  this  noise  is  troublesome  to  the  com¬ 
munity.  Fortunately,  the  characteristics  of  the  engine  as  a  source  of  community 
noise  can  be  fairly  easily  summarized.  They  have  to  do  with  overall  power, 
directivity,  and  spectra  as  these  are  related  to  engine  power  and  geometry  [16] . 
Unhappily,  the  loads  that  the  engine  noise  places  on  the  vehicle  are  not  so  simply 
related  to  these  parameters,  but  depend  on  details  of  the  vehicle  and  engine  con¬ 
figuration  in  a  fairly  complicated  way. 

The  long-range  propagation  of  engine  noise  has  been  reasonably  well  docu¬ 
mented,  and  will  not  be  pursued  here.  There  has  also  been  a  fair  amount  of  work 
on  the  structural  loads  induced  by  subsonic  jets;  some  of  it  is  documented  in  the 
open  literature,  but  much  of  it  resides  in  the  data  files  of  airframe  manufacturers. 
The  structural  loads  are  reviewed  in  the  following  paragraphs. 

The  determination  of  acoustic  loads  for  airframe  applications  is  perhaps  best 
exemplified  by  the  work  of  Clarkson.  In  two  particular  studies,  on  the  geometric 
zone  pressure  field  of  the  Comet  engine  [17,  p.  22]  and  in  a  mockup  of  the 
Caravelle  tail  section  [17,  p.  27],  Clarkson  has  obtained  data  on  the  scaling 
laws  for  pressure  correlations  between  jets  of  varying  sizes.  Using  these  correla¬ 
tions,  he  has  also  been  able  to  determine  the  boundary  between  the  nonradiated 
hydrodynamic  field  of  the  jet  and  the  geometric  field. 

The  induction  (hydrodynamic)  field  of  the  subsonic  jet  extends  about  2 
diameters  outward  from  the  jet  boundary.  In  this  region,  correlation  patterns  of 
pressure  are  similar  to  those  of  boundary  layer  pressure  fluctuation,  except  that 
the  dependence  of  convection  velocity  on  frequency  is  much  greater,  the  de¬ 
pendence  being  approximately 


Uc~f'n- 

This  is  a  much  sponger  relationship  than  similar  ones  derived  for  a  turbulent 
boundary  layer  [18] . 

The  intensity  and  spectra  of  pressure  fluctuations  in  the  induction  and  geo¬ 
metric  fields  have  not  been  fully  studied.  Most  workers  assume  that  the  induc¬ 
tion  field  is  proportional  to  the  dynamic  pressure  of  the  jet, 
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and  that  the  rms  pressure  in  the  geometric  field  [19]  (which  extends  outward 
about  10  diameters  from  the  jet  boundary  and  20  diameters  downstream)  varies 
with  jet  parameters  in  the  same  way  that  the  far-field  radiation  does.  These 
regions  are  summarized  in  Fig.  1. 

In  the  subsonic  jet,  the  total  radiated  power  rises  as  a  very  high  power  of  the 
jet  efflux  velocity  U.  If  the  jet  is  relatively  clean  aercdynf.mically,  the  law  is 
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Fig.  J .  Sketch  of  zone  geometry  for  subsonic  jet,  showing  induction, 
Fresnel  and  fraunhofer  radiation  regions  [  17]  .* 


Urad-6  X  lCr*pD*U*lcs.  (2.1) 

where  p  is  the  density  of  gas  in  the  jet,  D  is  the  nozzle  diameter,  and  c  is  the  speed 
of  sound  in  the  surrounding  medium. 

The  spectrum  of  overall  radiated  power  is  shown  in  Fig.  2  [16] .  The  spectrum 
can  generally  be  normalized  by  a  Strouhal  frequency  U/D.  Thus,  the  geometric 
field  magnitude  and  spectrum  can  usually  be  adjusted  to  predict  the  sound  from 
jets  having  differing  efflux  velocities  and  sizes,  providing  that  they  are  subsonic, 
“clean,”  and  that  one  is  comparing  spectra  at  geometrically  similar  observation 
positions. 

Supersonic  Jets 

When  the  flow  issuing  from  the  jet  nozzle  is  faster  than  the  speed  of  sound  in 
the  surrounding  medium,  the  efficiency  of  the  turbulence  sound  sources  will  be 
increased.  If  the  jet  is  hot,  as  it  usually  is  in  cases  of  interest  for  jet-poweied 
vehicles,  the  ilow  can  be  subsonic  with  respect  to  the  jet  itself,  but  supersonic 
with  respect  to  the  ambient  air.  In  this  case,  the  flow  in  the  jet  will  decelerate 
smoothly.  When  the  (low  in  the  jet  is  also  supersonic  with  respect  to  itself,  then 
strong  shocks  may  occur  as  it  decelerates. 

♦Figure  1  in  “Scaling  of  the  Near  Field  Pressure  Correlation  Patterns  around  a  Jet  Exhaust," 
by  B.  L.  Clarkson,  in  Acoustical  Fatigue  in  Aerospace  Structures,  edited  by  Walter  J.  Trapp 
and  Donald  M.  Forney,  Jr.  Copyright  ©  1965  by  Syracuse  University  Press,  Syracuse,  New 
York.  By  permission  of  the  publisher. 
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DIMENSIONLESS  FREQUENCY 


Fig.  2.  Normalized  power  spectrum  for  subsonic  jets  (16). 


Figure  3  is  a  Schlieren  photograph  of  a  supersonic  jet  showing  a  complicated 
pattern  of  sound  emission  from  the  supersonic  region  as  Mach  wavelets,  from  the 
shock  region,  and  from  the  subsonic  region.  The  relative  strength  of  each  portion 
of  the  jet  as  a  source  region  remains  to  be  clarified,  as  does  the  nature  of  the 
sources.  Nevertheless,  a  fair  amount  of  data  has  been  collected  from  jet  and 
rocket  engine  firings,  and  i‘  will  be  useful  to  review  some  of  these  and  how 
they  arc  applied  to  define  environmental  levels. 

Early  studies  of  model  supersonic  jets  were  carried  out  for  the  purposes  of 
understanding  their  noise  production-how  the  frequency  spectra  of  radiated 
power  scale  with  efflux  velocity  and  jet  diameter,  how  the  overall  radiated  power 
relates  to  jet  power,  and  how  the  distribution  of  radiated  power  depends  on  the 
angle  from  the  jet  axis.  Kor  these  purposes,  the  measurement  of  noise  pressures 
in  the  neighborhood  of  '.lie  jet-whether  they  be  the  nonradiated  induction  field 
or  the  radiated  but  complex  geometric  field  (or  Fresnel  zonc)-rcprescnts  an 
unnecessary  complication.  Thus,  much  of  the  information  on  rocket  noise  is  of 
little  value  in  estimating  the  loads  on  the  launch  vehicle  or  other  nearby  structures. 

We  can  see  examples  of  this  by  comparing  estimates  of  rocket  noise  on  a 
vehicle  with  published  data  on  flight  vehicles.  An  early  estimate  of  vehicle  en¬ 
vironment  was  given  by  Bics  and  Frankcn  [20] .  They  also  extrapolated  far-field 
measurements  to  the  vehicle  locations.  In  general,  these  estimates  predict  a 
frequency  maximum  in  the  excitation  that  is  too  low.  They  also  tend  to  be 
rather  low  in  amplitude. 

A  revision  to  the  early  estimates  was  made  in  1963  by  Frankcn  and  Wiener 
[21].  This  procedure  was  based  on  measurements  of  acoustic  spectra  on  launch 
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vehicles-Jupiter,  Atlas,  Titan,  and  Saturn  I  static  firings.  More  recent  data  on 
Saturn  launch  vehicles  tend  to  confirm  the  general  features  of  the  latter  esti¬ 
mates,  except  that  (a)  the  levels  do  not  diminish  as  quickly  as  anticipated  as  one 
proceeds  from  the  booster  to  the  payload  sections  and  (b)  the  spectrum  does 
not  undergo,  in  the  same  progression,  as  large  a  shift  to  lower  frequencies  as 
anticipated. 

The  conclusion  that  one  is  forced  to  draw  from  the  foregoing  is  that  reasonable 
estimates  of  vehicle  loads  due  to  acoustic  noise  from  its  jets,  subsonic  or  super¬ 
sonic,  must  at  present  be  based  on  direct  measurements  of  the  loads  on  similar 
or  scaled  vehicles.  Theoretical  notions  of  sound  generation  and  propagation  are 
useful  as  guides  to  scaling  the  frequency  scale  and  overall  level,  but  they  are  of 
little  value  at  present  in  predicting  the  levels  on  a  new  vehicle  configuration. 

This  rather  empirical  state  of  affairs  in  acoustic  load  estimation  is  repeated 
when  one  considers  the  other  flight  loads-TBL,  separated  flows,  wakes,  base 
pressure  fluctuations,  oscillating  shocks,  etc.  Paradoxically,  however,  the  extra¬ 
polation  of  laboratory  and  field  data  to  new  situations  for  these  loads  is  more 
straightforward  than  for  the  sound  field,  although  this  may  be  an  illusion  main¬ 
tained  by  the  paucity  of  published  data  on  aerodynamic  loads. 


NASA  Langley  Research  Center  photograph. 
Fig.  3.  Supersonic  jet  showing  Mach  wave  formation. 
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2.2  Inertial  Aerodynamic  Loads 

The  aerodynamic  loads  that  we  shall  discuss  are  due  essentially  to  the  forward 
motion  of  the  vehicle  and  may  be  broadly  categorized  into  two  classes.  The  first 
we  term  inertially  induced;  these  tend  to  be  proportional  to  the  dynamic  head  of 
the  flow  about  the  vehicle.  Examples  are  the  TBL,  separated  and  wake  flows  and 
base  pressure  fluctuations.  The  second  class  of  loads  is  compressibility  induced. 
Examples  of  these  are  transonic  and  supersonic  oscillating  shocks  and  cavity 
resonances. 

Based  on  wind  tunnel  observ  lions  and  flight  data,  the  locations  and  magni¬ 
tudes  of  the  major  fluctuating  aerodynamic  environments  are  fairly  well  known. 
A  summary  of  these  environments  and  their  occurrence  as  a  function  of  geom¬ 
etry  and  speed  is  shown  in  Fig.  4. 

An  attached  TBL  is  a  major  environment  over  large  parts  of  the  vehicle  for 
a  substantial  portion  of  the  atmospheric  flight  profile.  This  is  true  for  both  air¬ 
craft  and  spacecraft.  The  prediction  of  the  TBL  load  is  largely  empirical  and  is 
based  on  wind  tunnel  data  and  flight  data  for  both  aircraft  and  spacecraft.  A 
recent  report  by  Bies  [22]  summarizes  much  of  the  available  TBL  pressure  data 
for  aircraft,  spacecraft,  and  wind  tunnels.  In  a  subsonic  wind  tunnel,  the  overall 
fluctuating  pressure  coefficient  Cpnm  is  nearly 

Cpn„,  =  0.006.  (2.2) 

As  the  Mach  number  increases,  this  coefficient  tends  to  decrease,  and  for 
Reynolds  numbers  greater  than  107  (based  on  distance  from  the  leading  edge), 
the  reduction  is  nearly  10  dB  by  the  time  the  Mach  number  has  increased  to  3 
or  4. 

The  fluctuation  pressure  coefficients  for  aircraft  in  regions  of  smooth  flow 
appear  to  be  equal  to  those  for  flai  plates  in  wind  tunnels.  For  spacecraft,  there 
are  reasons  where  the  flow  is  disturbed,  not  separated.  In  such  regions  we 
may  therefore  prefer  to  use  the  coefficient 

Cprms  =  0.02.  (2.3) 

Additional  flight  data  will  improve  estimates  of  these  parameters,  but  these  will 
serve  satisfactorily  for  most  purposes. 

The  overall  pressure  fluctuation  is  of  interest  but  it  must  be  supplemented  by 
a  knowledge  of  the  spatial  and  temporal  corielations  (or  spectra)  in  order  to 
predict  response.  To  convert  from  the  fixed-position  single  microphone  measure¬ 
ment  to  a  determination  of  the  correlation,  a  model  of  the  pressure  correlations 
must  be  assumed.  A  particular  model  of  the  correlation  that  has  enjoyed  wide 
usage  is  the  following  [23] : 
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Fig.  4.  Fluctuating  aerodynamic  loads  depending  on  vehicle  location  and  Mach 
number  [2], 
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(p(x.  t)p(x,  /')>  =  <p2>0p($,,  £3,  o 

=  <^)0i(^i-^r)03a3)0m(r)  (2.4) 

where  £,  =  xt  —  x\,  £3  =  x3  —  .v'3,  t  =  t  —  t'.  Equation  (2.4)  represents  a 
correlation  field,  convected  in  the  .v,  (downstream)  direction  at  a  speed  Uc  and 
changing  in  time  so  that  the  temporal  correlation  in  this  moving  frame  is 
0m(r).  The  change  in  correlation  as  one  traverses  across  the  flow  is  given  by  $y 
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The  form  of  Eq.  (2.4)  is  a  reflection  of  the  way  in  which  TBL  pressure  data 
are  analyzed.  A  measurement  of  the  autocorrelation  at  one  point  is  given  by 
setting  | !  =  |3  =  0: 

<t>pif)  =  <t>y{VcT)<}>m{r).  (2.5) 

It  is  usually  assumed,  on  the  basis  of  empirical  evidence,  that  <Pm(r)  is  much  more 
slowly  varying  than  is  <pv  Thus,  at  least  for  small  values  of  £j  (or  Uct),  the 
fixed  microphone  measurement  <j>p  gives  a  good  indication  of  the  form  of  <j>1 .  It 
is  usually  a  function  that  decays  rapidly  in  «pace  and  becomes  negative  for  large 
values  of  the  argument.  There  has  been  a  long  standing  discussion  among  fluid 
dynamicists  as  to  whether  or  not  the  integral  /^t(f)^  has  a  non-zero  value,  but 
it  is  a  point  of  relatively  small  importance  for  structural  response  [24] . 

The  transverse  spatial  correlation  03(£3)  may  be  measured  by  setting  r  =  0 
and  -  0,  and  measuring  the  spatial  crosscorrelation  transverse  to  the  direction 
of  flow.  In  this  case,  there  is  fairly  good  correlation  of  results  between  experi¬ 
menters,  all  of  them  more  or  less  agreeing  on  the  form  [5;  23,  p.  15]: 

^3(^3)  ~  exP  [~ £3/^3 1  •  (2.6) 

There  are  also  data  taken  on  the  more  general  correlation  $p(£, ,  £3 , 0).  It  does 
not  completely  support  the  factoring  of  the  spatial  correlation  into  separate 
downstream  and  cross-stream  components.  Almost  all  applications  of  loads  in¬ 
formation  at  the  present  time  do  assume,  however,  that  this  factorization  can 
be  carried  out  [25] . 

The  experimental  determination  of  0,„(t)  is  suggested  by  Eq.  (2.4).  Since 
0j(O)  =  <A3(0)  =  I  ,<t>m  can  be  found  from 

0m(r)  =  fytt,  -t/c^3=0,r).  (2.7) 

The  experimental  process  consists  of  generating  a  series  of  curves  <Pp(%lt  0,  r), 
where  £,  is  a  fixed  parameter  and  r  is  the  abscissa.  The  function  0,„ (r)  is  the 
envelope  of  these  curves,  and  tire  convection  speed  is  the  ratio  of  the  £  value  for 
any  particular  curve  to  the  time  delay  r  at  the  tangent  point  [26] . 

The  preceding  discussion  is  oversimplified,  as  anyone  who  has  made  such 
experiments  can  testify.  In  fact,  the  values  of  convection  speed  as  obtained  for 
the  different  values  of  £,  will  vary  and  the  shape  of  <f>p  will  change  as  a  function 
of  r  in  a  way  that  is  not  totally  explainable  from  the  shape  of  the  modulation 
<Pm(r).  Much  time  and  effort  have  gone  into  trying  to  get  better  correlation 
models,  and  with  some  success.  Let  us  note,  however,  that  environmental  esti¬ 
mation  is  a  systems  problem.  It  has  not  been  demonstrated  that  many  of  the 
finer  points  in  the  correlation  determinations  have  any  significant  effect  on 
structural  response,  and  in  some  cases,  one  can  demonstrate  that  there  is  very 
little  effect  [23,  Fig.  1] .  An  example  is  the  insensitivity  of  response  of  a  high¬ 
speed  vehicle  to  the  moving  axis  temporal  correlation. 
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The  estimation  of  TBL  loads  relies  heavily  on  wind  tunnel  measurements,  but 
increasingly,  flight  data  are  also  playing  a  part  in  the  loads  estimates.  By  far,  the 
bulk  of  the  data  is  of  single-microphone  temporal  autocorrelations  or  of  its 
frequency  transform,  the  power  spectral  density  (PSD).  In  a  few  cases,  spatial 
crosscorrelations  have  been  measured,  and  in  general,  they  support  the  conclusions 
drawn  from  the  wind  tunnel  measurements,  except  for  the  increase  in  overall 
level  as  noted  above.  Flight  data  have  been  obtained  on  a  glider  wing  [27] ,  on 
commercial  [28]  and  military  aircraft  [29] ,  and  for  spaceciaft  [30] ,  launch 
vehicles  [31]  and  experimental  military  payloads  [32].  The  TBL  environment 
can  be  summarized  as  follows: 

1.  The  overall  fluctuating  pressure  varies  with  dynamic  head  and  Mach  num¬ 
ber,  from  approximately  0.02<2  at  M  =  1  to  0.006  at  M  =  4  for  spacecraft  and 
0.006  for  aircraft.  The  mission  profile  will  determine  the  coefficient  appro¬ 
priate  to  the  period  of  high  TBL  excitation. 

2.  The  pressure  disturbances  are  convected  along  the  vehicle  surface  at  60  to 
80  percent  of  the  local  free-stream  flow  speed. 

3.  The  spatial  scale  of  the  disturbance  in  both  the  down-  and  cross-stream 
directions  is  approximately  6*,  where  5*  is  the  displacement  thickness  of  the 
boundary  layer.  The  actual  spectra  of  0,  and  <p3  are  required  for  response 
estimations,  and  these  have  the  fairly  universal  forms  shown  in  Figs.  5  and  6. 

Separated  Flows 

Separated  flow  fields,  when  they  occur  on  a  flexible  structural  surface,  can 
produce  a  significant  part  of  the  response  due  to  the  fluctuating  pressure  en¬ 
vironment.  In  this  category  we  can  place  wakes  from  protuberances  and  other 
nonfaired  structural  elements  on  the  forward  parts  of  the  vehicle,  shock-induced 
separation  on  transonic  and  supersonic  vehicles  near  places  where  there  is  a  rapid 
rate  of  change  of  vehicle  cross-sectional  area,  and  stalled  flow  regions  such  as  at 
base  separation  or  on  separation  airfoils  at  large  angles  of  attack. 

In  the  following  chapter,  we  discut  s  methods  for  using  loads  information  to 
estimate  structural  response.  Suffice  it  to  say  here  that,  in  general,  the  ratio  of 
response  to  excitation  in  any  frequency  band  depends  on  the  type  of  excitation. 
A  particular  pressure  PSD  measured  at  a  position  on  the  structure  will  produce 
differing  response  PSD’s  if  the  pressure  is  due  to  a  sound  field  or  if  it  arises  from 
a  TBL.  We  must,  therefore,  know  more  than  the  PSD-we  must  know  what  kind 
of  flow  induces  that  pressure  field  before  we  can  calculate  response.  Unfortu¬ 
nately,  only  a  very  few  correlation  studies  have  been  made  for  separated  flows. 
A  greater  degree  of  guesswork  (the  euphemism  is  extrapolation)  is  required  in  the 
response  computation  than  when  the  excitation  arises  from  an  acoustic  field 
or  a  TBL. 

A  useful  and  simple  assumption  is  that  a  region  of  separated  or  wake  flow  is 
like  a  very  thick  TBL,  with  a  convection  speed  approximately  50  to  60  percent 
of  the  free  stream  (or  vehicle)  speed  [2].  The  intensity  of  the  pressure 
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fluctuations  in  separated  regions  on  the  forward  part  of  the  vehicle  associated 
with  transonic  and  supersonic  shocks  is  approximately  [2] 


Pnm  -  0.05 Q.  (2.8) 

In  regions  where  wake  flow  impinges  on  the  vehicle,  the  pressure  is  more 
nearly  [2] 


Pnns  —  0. 1 Q. 


(2.9) 


FLIGHT  VEHICLE  NOISE  ENVIRONMENT  1 5 

There  have  been  some  measurements  of  both  pressure  fluctuation  and  corre¬ 
lation  for  separated  flow  over  a  delta  wing  aircraft  at  Southampton  [33] .  In  this 
case,  the  overall  pressure  fluctuation  was  nearly 

Prms  ~~  0.02 Q.  (2.10) 

The  pressure  fluctuations  did  show  a  convection  effect,  but  the  convection  had  a 
large  component  transverse  to  the  main  direction  of  flow. 

For  purposes  of  response  estimation,  one  can  assume  that  the  convected  loads 
due  to  separated  flows  are  completely  analogous  to  a  TBL  with  a  larger  charac¬ 
teristic  dimension.  For  a  wake,  the  radius  dimension  of  the  obstacle  producing 


Fig.  6.  TBL  pressure  wave  number  spectrum  in  cross-stream  direction  1 23  J . 
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the  wake  may  take  the  place  of  the  displacement  thickness  8*.  In  supersonic 
or  transonic  shock-induced  separation,  it  is  a  multiple  of  the  displacemem  thick¬ 
ness  of  the  incoming  TBL-the  factor  is  usually  2  or  3  [2] . 

At  the  base  of  the  vehicle,  or  at  a  more  forward  part  of  the  vehicle  that  is 
strongly  reentrant,  strong  flow  separation  is  also  accompanied  by  pressure  fluc¬ 
tuations.  Such  loads  are  called  base  pressure  fluctuations  (BPF).  Unfortunately, 
the  only  data  existing  give  only  PSD  information  for  a  single  microphone  sensor 
[34] .  There  is  essentially  no  information  available,  therefore,  on  the  spatial 
correlation  or  the  presence  (or  absence)  of  convection  effects.  Also,  the  data  are 
for  a  subsonic  case.  These  data  do  suggest,  however,  two  important  features  of 
this  environment.  First,  the  overall  pressure  fluctuations  are  closer  in  magnitude 
to  the  attached  TBL  (Eq.  2.3)  rather  than  the  separated  flow  values  in  Eqs.  (2.8) 
through  (2.10).  Secondly,  the  PSD  was  essentially  the  same  shape  as  Fig.  5  for 
the  TBL,  except  that  the  diameter  of  the  vehicle  should  replace  the  displacement 
thickness  as  a  characteristic  dimension  and  co/£/  replaces  . 

The  question  of  convection  has  not  been  studied  experimentally,  but  it  seems 
logical  that  convection  should  be  absent  in  BPF  loads.  If  this  is  the  case,  spatial 
correlation  experiments  will  still  be  required  to  determine  the  form  of  ,|3 ,0). 
In  the  absence  of  such  data,  one  may  assume  a  correlation  form  for  estimation 
purposes.  The  replacement  of  the  environment  with  a  sound  field  is  one  example. 
Another  possibility  is  to  assume  a  relation  between  spatial  scale  and  correlation 
that  is  logical  and  convenient,  for  example,  the  form  of  Eq.  (2.6)  with  L  =  D, 
the  vehicle  (or  base)  diameter. 

Base  Pressure  Fluctuations 

Base  pressure  fluctuations  (BPF)  occur  in  regions  of  separated  flow  on  the  aft 
section  of  a  vehicle  as  shown  in  Fig.  7.  Pressure  fluctuations  in  this  region  of  the 
vehicle  may  be  expected  to  be  somewhat  different  from  those  at  more  forward 
positions  on  the  vehicle  for  two  reasons:  (a)  there  is  no  reason  to  expect  signif¬ 
icant  convection  effects  for  BPF,  and  (b)  the  spatial  scale  is  set  by  the  overall 
vehicle  diameter. 

There  are  only  two  pieces  of  data  on  BPF  that  are  generally  available.  Onv 
is  the  previously  mentioned  subsonic  data  taken  by  Eldred  in  a  wind  tunnel  [34] , 
which  by  showing  that  the  frequency  spectrum  of  BPF  is  similar  to  that  of  the 
TBL  in  the  attached  regions  of  flow,  except  that  the  vehicle  diameter  replaces  the 
displacement  thickness  as  the  characteristic  dimension,  generally  leads  to  a  spec¬ 
trum  with  much  greater  low-frequency  content.  Also,  since  the  overall  rms 
pressure  coefficient  in  the  base  is  very  nearly  equal  to  that  in  the  region  of 
attached  flow,  the  low-frequency  loading  is  very  much  greater  on  the  base  of  the 
vehicle  than  on  its  sides— the  opposite  is  true,  of  course,  for  high-frequency 
excitation. 

The  second  piece  of  data  is  that  inferred  by  Bcncdetti  and  Smith  [32]  from 
base  vibration  and  internal  microphone  data  on  an  experimental  reentry  vehicle. 
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Fig.  7.  Reentry  vehicle  aerodynamics  showing  regions  of  at¬ 
tached  (TBL)  and  separated  (BPF)  pressure  loads. 


Unfortunately,  the  amount  of  data  reported  do  not  seem  sufficiently  conclusive 
to  allow  comparison  with  the  subsonic  data  of  Eldred. 

Very  little  is  known  about  the  spatial  correlation  of  BPF.  An  appealing 
assumption  is  to  use  the  transverse  spatial  correlation  of  the  TBL,  replacing  the 
displacement  thickness  5*  by  A  the  vehicle  diameter: 

<p  ~  exp  ( —  r/D). 

There  is  no  experimental  support  for  this,  but  it  can  be  useful  in  estimating 
response  of  base  structures  and  probably  produces  a  result  in  the  .ight  ballpark. 

2.3  Compressible  Aerodynamic  Loads 

Tire  most  important  noninertial  aerodynamic  load  on  a  flight  vehicle  arises 
from  the  existence  of  shock  waves.  When  a  shock  occurs,  it  is  rarely  absolutely 
steady,  and  the  fluctuation  in  its  position  will  cause  pressure  variations  on  the 
vehicle  surface,  generally  at  relatively  low  frequencies.  The  origin  of  the  fluc¬ 
tuations  is  not  clear-but  they  arc  thought  to  be  associated  with  the  upstream 
turbulent  boundary  layer  and  region  of  separated  flow  downstream.  Figure  8 
is  a  photograph  of  a  space  vehicle  model  in  a  wind  tunnel  that  demonstrates  rather 
clearly  the  incoming  TBL,  the  shock,  and  the  region  of  separated  flow  [35]. 

Measurements  of  fluctuating  pressures  in  the  region  of  an  attached  shock 
generally  show  augmented  low-frequency  pressure  levels.  The  oscillation  of  the 
shock  is  one  source  of  low-frequency  loads,  but  the  interaction  of  the  shock  with 
the  inflowing  TBL  will  cause  acoustic  pressure  fluctuations,  both  radiated  and 
nonradiated.  Lowson  has  shown  theoretically  that  quite  large  pressure  fluctua¬ 
tions  can  be  explained  from  this  mechanism  [36] .  It  is  not  clear  at  this  time, 
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however,  whether  or  not  Lowson’s 
model  can  account  for  the  apparent 
magnitude  of  the  random  displace¬ 
ment  of  the  shock  front  or  the  ot- 
served  spectrum.  In  any  case,  this 
part  of  the  TBL-shock  interaction 
must  be  present  and  will  represent  a 
contribution  to  the  loads. 

Low-frequency  fluctuations  will 
also  arise  from  the  shock-separated 
flow  interaction.  The  higher  pres¬ 
sures  in  the  downstream  side  of  the 
shock  cause  flow  separation.  The 
separated  flow  region  behaves  like  a 
TBL  with  a  very  thick  displacement 
thickness.  This  separated  flow  gen¬ 
erates  low-frequency  pressures.  Also, 
the  high  downstream  pressure  is  able 
to  penetrate  the  boundary  layer  and 
cause  it  to  thicken  upstream.  This 
thickening  reduces  the  shock  strength 
and  causes  the  TBL  to  thicken  even 
more-resulting  in  the  displacement 
of  the  shock  to  a  new  stable  position 
upstream  [37] .  Once  the  shock  moves 
to  this  new  position,  however,  the 
downstream  conditions  cannot  main¬ 
tain  it  there,  and  it  returns  whence  it 
came.  And  so  it  goes,  back  and  forth. 

This  latter  description  seems  to  be 
drawn  from  analyses  of  certain  buf¬ 
feting  problems  [38].  It  may  be 
adequate,  however,  to  call  out  the 
two  phenomena,  the  shock  motion 
and  the  separated  flow,  as  significant  loads.  The  separated  flow  has  already 
been  treated;  here  we  will  treat  the  shock  motion. 

Chandiramani  ct  al  [2]  describe  loads  due  to  supersonic  shock  oscillation,  of 
the  kind  that  occur  ahead  of  a  flare  in  the  structure,  by  the  displacement  of  the 
shock  y(t)  and  the  pressure  differential  across  the  shock  Ap.  The  displacement 
is  assumed  to  be  a  narrow-band  process  centered  about  frequency  fox  with 
bandwidth  foscl 4.  The  magnitude  Ap  is  the  pressure  differential  obtained  from 
the  supersonic  aerodynamics  of  the  vehicle.  The  rms  shock  displacement  o  is 
taken  to  be  a  fraction  1/8  of  the  upstream  distance  L  from  the  flare  to  the  shock. 
The  frequency  fosc  is  obtained  from 


Fig.  8.  Space  vehicle  model  in  supersonic 
flow  showing  shock-induced  separation  for¬ 
ward  of  flare  [35]. 
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fosc~  0.06  U/L, 

where  U is  the  speed  of  the  vehicle.  The  spectrum  of  the  signal  that  this  produces 
in  a  microphone  is  shown  in  Fig.  9.  Note,  however,  that  this.is  not  the  same  as 
the  effective  spectrum  of  generalized  force  on  any  particular  mode  of  the  struc¬ 
ture.  We  shall  take  up  the  response  problem  in  the  following  chapter.  A  sum¬ 
mary  of  the  estimation  procedure  for  oscillating  shock  pressure  spectra,  their 
period  of  occurrence  during  the  flight,  and  their  location  on  the  vehicle  is  by 
by  Chandiramani  et  al.  [2] . 

The  preceding  discussions  give  as  complete  a  picture  as  possible  for  the  cur¬ 
rent  status  of  flight-induced  fluctuating  pressure  environments.  More  details  of 
these  environments,  particularly  as  they  occur  on  space  vehicles,  can  be  found  in 
Ref.  2.  Nevertheless,  it  is  clear  that  there  are  many  gaps  in  our  knowledge  of 
flight-induced  loads.  It  is  important,  however,  that  we  not  take  data  for  its  own 
sake,  but  that  its  intended  use  be  considered  in  the  experimental  design. 

Far  too  often  great  stress  has  been  placed  in  accumulating  data  that  have 
turned  out  to  be  peripheral  or  beside  the  point.  A  good  example  of  this  is  the 
concern  whether  the  correlation  0,  (^, )  in  a  TBL  does  or  does  not  have  a  finite 
correlation  area.  Structural  analysts  would  be  wise  to  give  their  requirements  to 
the  aerodynamicists  early  in  the  loads  accumulation  program.  In  the  following 
chapter,  we  have  a  look  at  response  prediction  models.  These  will  indicate  the 
loads  data  needed  to  predict  response. 


CHAPTER  3 

ESTIMATION  OF  VEHICLE  RESPONSE 


3.1  Introduction 

The  design  engineer  and  the  environmental  or  test  engineers  are  interested  in 
the  response  of  a  vehicle  to  in-flight  loads  for  slightly  differing  reasons.  The  test 
engineer  is  primarily  concerned  with  the  fact  of  vibration-Acnv  much  vibration, 
when,  and  where.  The  design  engineer  wants  the  how  much  information,  so  that 
he  can  assess  the  integrity  of  the  vehicle.  But  he  is  also  concerned  with  the  why 
of  the  response,  since  he  wants  to  know  what  he  can  change  in  the  event  that  the 
vibration  is  excessive,  or  to  anticipate  the  effects  of  structural  modifications. 

The  procedures  for  response  estimation  reflect  this  dual  purpose  to  a  degree, 
but  it  is  safe  to  say  that  the  test  engineer  has  been  better  served.  Fairly  empir¬ 
ical  estimation  procedures  are  in  common  use,  but  there  has  been  an  effort  to 
develop  methods  that  relate  the  response  more  directly  to  structural  parameters. 
Some  of  these  methods  still  rely  on  empiricism  for  portions  of  the  prediction 
procedure-we  shall  term  these  semi-empirical.  Quite  recently,  some  procedures 
have  been  put  forward  that  rely  much  less  on  empiricism  and  much  more  on 
analysis  of  structural-environmental  interactions.  We  refer  to  these  as  rational 
methods. 

The  goals  of  response  estimation  are  (a)  the  providing  of  information  on  the 
vibration  of  the  external  structure,  >ts  frequency  and  spatial  distribution  and  how 
these  change  with  time;  and  (b)  the  transmission  of  this  vibrational  energy  to  in¬ 
ternal  structures  and  acoustic  spaces.  As  the  various  procedures  are  described, 
the  extent  to  which  each  meets  these  goah  will  be  discussed. 

3.2  Empirical  and  Semi-empirical  Methods 

The  most  direct  approach  to  response  estimation  is  to  measure  the  vibration 
and/or  sound  levels  experienced  by  a  vehicle  during  its  flight,  and  then  to  use 
these  levels  to  predict  the  vibration  of  a  panel,  or  black  box,  that  is  located  on  or 
.thin  a  similar  vehicle  flying  a  similar  trajectory.  Even  such  a  straightforward 
procedure  has  difficulties  that  are  worth  mentioning. 

A  major  difficulty  derives  from  the  method  of  data  collection.  Space  vehicle 
systems  thus  far  use  telemetry  rather  than  on-board  recorders  to  gather  most  of 
the  data.  In  many  cases  the  vehicle  is  not  recovered;  in  others  the  weight  of  a 
recorder  isjudged  too  great  for  inclusion  in  the  payload.  The  frequency  response 
of  a  telemetry  channel  is  limited- a  bandwidth  up  to  a  few  kilohertz  at  the  most- 
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and  the  number  of  channels  is  limited.  Thus,  the  relayed  information  does  not 
usually  contain  much  high-frequency  data. 

Recorders  can  be  carried  in  aircraft,  and  the  frequency  range  of  available  data  is 
greater  than  it  is  for  spacecraft.  A  difficulty  shaied  by  both,  however,  is  that  due  to 
the  size  of  these  vehicles  and  the  many  kinds  of  structures  involved,  the  number  of 
observation  locations  is  not  adequate  for  good  determination  of  response  or  for  ade¬ 
quate  definition  of  the  subsystem  environment.  Nevertheless,  with  repeated  flights, 
a  good  library  of  flight  data  car.  be  built  up,  and  the  vibration  or  internal  sound  lev¬ 
els  at  a  great  number  of  locations  during  a  flight  can  be  predicted  for  the  vehicle. 

The  methodology  involved  in  this  data  collection  is  rather  straightforward.  The 
process  of  converting  the  data  into  a  prescribed  vibration  environment  may  be  rather 
more  involved,  however.  Typically,  the  output  of  a  three-axis  accelerometer  located 
at  some  sensitive  position  is  monitored  as  a  function  of  time.  This  signal  is  ana¬ 
lyzed  in  frequency  bands,  and  the  spectrum  at  different  times  is  plotted.  The  levels 
are  likely  to  be  large  during  specific  portions  of  the  flight  (launch,  transonic,  or 
max  Q).  The  spectra  for  these  portions  of  the  flight  are  then  converted  into  test 
power  spectral  densities,  usually  prescribed  by  the  slope  and  magnitude  of  the  spec¬ 
trum  over  frequency  intervals.  It  is  in  the  conversion  of  the  environmental  data  to 
test  spectra  that  most  of  the  vagaries  arise. 

A  typical  example  of  reduced  vibration  data  is  shown  in  Fig.  10.  Such  a  spec¬ 
trum  is  too  complicated  to  prescribe  as  a  test  spectrum,  and  in  addition,  it  is  only 
a  typical  one-other  flights  and  other  locations  will  produce  similar,  but  slightly 
differing  spectra  (variations  in  band  levels  of  the  order  of  10  dB  are  common).  The 
simplest  procedure  is  to  envelope  this  spectrum  with  a  spectrum  that  can  be  speci¬ 
fied.  Curve  A  in  Fig.  10  is  an  example.  In  most  cases  this  spectrum  will  produce  a 
conservative  test;  i.c.,  the  test  is  more  severe  than  the  environment.  As  an  alter¬ 
native,  curve  A  may  be  lowered  until  the  overall  acceleration  level  of  the  test  is  the 
same  as  that  of  the  environment,  and  the  new  test  spectrum  is  curve  B.  In  this  case, 
some  of  the  peaks  in  the  realized  environmental  spectrum  will  exceed  the  test 
spectrum,  and  the  engineer  may  be  concerned  that  a  resonance  in  the  tested  item 
at  such  frequencies  may  be  underexposed. 

In  addition  to  these  problems,  the  adequacy  of  the  environmental  data  may  be 
in  question.  If  several  flights  are  available,  a  composite  of  the  data  in  Fig.  10  can 
be  drawn  which  exceeds  80,  90,  or  1 00  percent  of  the  data.  This  in  turn  can  be 
converted  to  a  test  spectrum  by  an  enveloping  procedure. 

The  decision  process  by  which  one  converts  environmental  data  to  a  test  spec¬ 
trum  is  referred  to  as  test  philosophy.  To  the  cynic,  uninitiated  into  the  rites  of 
such  activities,  a  more  appropriate  word  might  be  theology.  In  many  cases,  the 
exact  procedure  is  unimportant-the  item  is  developed  and  tested  for  levels  above 
its  environment.  As  we  shall  see  in  the  following  chapter,  however,  this  is  far  from 
being  the  only  set  of  vagaries  in  the  system.  Despite  its  importance,  the  problem  of 
setting  test  levels  has  received  rather  little  research  effort.  It  could  well  benefit  from 
some  basic  studies  of  the  effect  of  fine-grained  spectral  variations  on  structural  or 
component  response. 
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Fig.  10.  Narrow-band  vibration  spectrum 
and  derived  (smooth)  test  specifications. 


The  Method  of  Mahaffey  and  Smith 

The  trajectory  (or  mission)  profile  of  a  flight  vehicle  tends  to  be  controlled  by 
its  design  and  configuration.  Some  variations  will  be  encountered.  The  trajectory 
of  a  launch  vehicle  will  depend  on  whether  the  flight  is  suborbital,  orbital,  or 
interplanetary.  An  aircraft  has  a  different  trajectory  in  landing,  takeoff,  recon¬ 
naissance,  or  attack  [39] . 

In  addition  to  trajectory  variations,  the  vehicle  itself  may  have  structural 
changes.  Its  engines  may  be  relocated,  it  may  be  shortened,  or  certain  sections 
may  have  smaller  geometrical  changes.  To  capitalize  on  information  previously 
gathered  on  the  vehicle,  several  methods  have  been  devised  to  allow  one  to  trans¬ 
late  the  old  data  to  new  circumstances.  One  of  these  is  the  so-called  Mahaffey- 
Smith  (M-S)  method.  We  will  illustrate  this  method  by  employing  the  same 
example  used  by  Mahaffey  and  Smith. 

The  initial  step  is  a  fairly  complete  sound  and  vibration  survey-in  this  in¬ 
stance  a  measurement  of  acceleration  and  sound  pressures  analyzed  in  octave 
bands  at  various  locations  on  the  B-58  airplane  shown  in  Fig.  1 1.  The  authors 
expressed  the  response  as  peak  acceleration,  which  is  three  times  the  rms  value 
the  response;  defining  the  peak  as  that  level  that  is  crossed  1  percent  of  the  time 
(which  is  consistent  with  the  response  being  Gaussian).  In  Fig.  12,  we  have  re¬ 
produced  a  set  of  the  data  for  the  600  to  1200  Hz  octave  band,  using  the  rms 
as  a  measure  of  the  acceleration. 

The  data  generally  appear  to  fit  a  linear  regression  line  rather  well  in  all  the 
octave  bands  studied.  Assuming  a  log-normal  distribution  for  the  rms  accelera¬ 
tions  for  a  given  sound  pressure  level  Lp  (vibration  data  in  complex  structures 
seem  to  fit  this  distribution  surprisingly  well  in  a  wide  variety  of  circumstances), 
the  variance  of  the  measured  data  is  noted  and  confidence  levels  for  each  band 
are  also  drawn,  as  shown.  In  general,  the  iorm  of  the  regression  line  is 
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Fig.  1 1.  Contours  of  constantly  on  B-58  airplane  for  600  to  1200 
Hz  octave  band  [40]  - 


La  =  MLp  +  A,  (3.1) 

where  A  is  a  measure  of  the  acceptance  of  the  structure  and  M  is  the  slope. 
Both  M  and  A  vary  from  one  frequency  band  to  the  next,  but  the  range  o CM  is 
from  0.5  at  the  lower  bands  to  0.7  in  the  upper  bands. 

The  observation  that  M  r  1  has  led  some  observers  to  suspect  nonlinear 
behavior  of  the  structure.  There  are  other  more  likely  reasons  for  a  lack  of 
proportionality  between  pressure  and  acceleration,  however.  One  of  these  is  the 
tendency  of  vibrational  energy  to  propagate  through  the  structure  from  a  region 
of  high  excitation  to  a  region  of  low  excitation.  Thus,  regions  of  high  excitation 
have  less  response  than  if  they  were  not  connected  to  oilier  structure,  while 
regions  of  lower  excitation  have  more.  Whether  mechanical  transmission  can 
account  for  the  behavior  of  the  field  data  has  not  been  determined  at  present. 

Finally,  one  derives  vibration  estimates  from  the  set  of  curves  such  as  those  in 
Fig.  1 2  by  measuring  the  sound  pressure  levels  Lp  on  a  section  of  the  vehicle  not 
previously  measured.  The  values  of  Lp  for  each  frequency  band  are  plotted,  and 
the  curves  like  Fig.  12  are  consulted.  The  appropriate  confidence  coefficient  is 
then  chosen,  and  the  corresponding  acceleration  level  in  each  band  is  read.  This 
set  of  levels  vs  frequency  becomes  the  response  estimate,  and  may  be  used  as  is, 
or  averaged  (smoothed)  in  some  way  to  serve  as  an  environmental  specification. 

It  is  relatively  easy  to  find  fault  with  the  M-S  procedure.  It  ignores  struc¬ 
tural  differences  throughout  the  vehicle,  it  is  not  specific  as  to  how  measure¬ 
ments  are  to  be  taken,  and  it  is  devoid  of  structural  or  geometric  scaling  param¬ 
eters.  For  these  reasons,  and  others,  it  has  on  occasion  been  misused.  It  is  more 
instructive,  however,  to  inquire  why  it  has  received  so  much  attention.  There 
are  at  least  three  reasons  why  this  is  so.  First,  the  M-S  procedure  is  relatively 
simple  to  carry  out.  It  requires  the  existence  of  a  similar  operational  vehicle,  and 
the  taking  of  a  fair  amount  of  sound  and  vibration  data  to  obtain  the  basic 
information.  The  gathering  of  this  information  is  straightforward,  however. 
Second,  it  generates  an  estimate  that  is  easy  to  interpret.  Finally,  the  estimate  of 
response  is  in  the  proper  form  for  application. 
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Franken  Estimate  for  Cylindrical  Vehicles 

The  M-S  procedure  does  not  take  into  account  the  effect  of  any  structural 
differences  on  the  response  at  two  different  portions  of  the  vehicle  exposed  to 
the  same  fluctuating  pressure  levels.  Franken  has  pointed  out  that  if  one  ex¬ 
amines  the  response  of  cylindrical  vehicles  to  acoustic  excitation,  then  the  ratio 
of  rms  acceleration  to  rms  pressure  reaches  a  maximum  near  the  ring  frequency 
fr  for  the  vehide,  and  has  a  rapid  falioff  above  fr  [41].  The  ring  frequency  is 
that  frequency  for  which  the  longitudinal  wavelength  in  the  skin  material  is  equal 
to  the  vehicle  circumference  2m: 

fr  =  Cillira.  (3.2) 

The  second  observation  is  that  the  response  appears  to  be  well  normalized  to 
the  limp  wail  response  of  the  vehicle,  the  response  that  the  pressures  would  in¬ 
duce  into  a  structure  of  the  same  surface  mass  density  but  having  no  rigidity. 
For  such  a  vehicle,  the  ratio  of  acceleration  to  pressure  is  constant 


Fig.  12.  Scatter  diagram  of  acceleration  levels  La  on  B-58  airplane 
plotted  against  pressure  levels  Lp  in  600  to  1200  Hz  octave  band  (401 . 
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WV/p2  =  1.  (3.3) 

Using  a  pressure  reference  of  0.0002  dyne/sq  cm,  the  acceleration  of  gravity  as  a 
reference  for  acceleration  level,  and  expressing  the  surface  density  of  the  panel 
W  in  pounds  per  square  foot,  one  gets  the  mass  law  response  formula 

La(ML)  +  20  log  W-LP  =  -  128.  (3.4) 

The  mass  law  acceleration  is  plotted  in  Fig.  13  along  with  the  cross-hatched 
band  representing  the  distribution  of  observed  response  on  Jupiter  and  Titan  I 
vehicles  as  analyzed  by  Franken.  We  see  that  the  actual  response  increases  from 
the  mass  law  base  line  to  approximately  8  dB  above  mass  law  at  the  ring  fre¬ 
quency  and  drops  off  to  mass  law  levels  at  a  few  octaves  above  fr. 


Fig.  13.  Normalized  acceptance  of  cylindrical  vehicles  of  diameter  D  and  surface  weight 
IV  (41], 


We  have  terme'1  Franken’s  method  semi-empirical  because  it  does  use  notions 
of  structural  detail  such  as  ring  frequency  and  surface  density.  It  nevertheless  is 
empirical  in  that  it  relies  exclusively  on  data  for  the  establishment  of  an  estimate. 
It  is  not  clear  at  this  point  whether  it  is  consistent  with  the  M-S  method.  It  may 
be,  since  in  an  aircraft,  regions  of  high  excitation  tend  to  have  smaller  diameter 
(see  Fig.  1 1)  and  the  change  in  fr  and  surface  density  as  one  proceeds  along  the 
fuselage  may  combine  to  give  results  such  as  those  in  Fig.  1 2.  The  available  data 
from  the  M-S  paper  do  not  permit  this  comparison. 
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Franken’s  procedure  has  been  very  useful  for  estimating  skin  vibration  for 
cylindrical  vehicles.  It  may  be  possible  to  develop  extensions  of  it  to  predict  the 
transmission  of  sound  and  vibration  to  interior  compartments  and  structures. 


Response  Estimation  Procedure  of  Eldred 

The  fact  that  Franken’s  empirical  estimate  is  8  dB  or  so  above  mass  law 
response  indicates  that  it  is  the  resonant,  damping  limited  response  of  the  struc¬ 
ture  that  is  dominant  near  the  ring  frequency.  In  a  series  of  reports  and  papers, 
Eldred  and  his  colleagues  have  attempted  to  use  this  observation  to  inject  addi¬ 
tional  structural  information  into  the  estimation  of  response  [42].  They  have 
done  this  by  developing  an  analogy  to  the  response  of  a  simple  resonator. 

Consider  a  resonator  with  a  resonance  frequency  f0,  mass  M  and  damping 
loss  factor  r\  where  t?  «  1.  Suppose  the  random  force  applied  to  the  resonator 
has  a  rather  smooth  frequency  spectrum  S(f).  Tire  mean  square  acceleration  of 
the  resonator  is  given  by  [43] 


S(fo)‘fo 

rjM2 


(3.5) 


In  analogy  with  this,  the  skin  acceleration  of  a  structure  of  area  A  and  surface 
density  IP  is  assumed  to  be  of  the  form 


Sp(f)f 
Vw2  ’ 


(3.6) 


where  Sp(f)  is  the  spectral  density  of  the  pressure  fluctuations.  The  parameter  j3 
is  a  measure  of  how  well  the  pressure  field  is  able  to  excite  the  niode(s)  that 
resonate  in  the  frequency  band  under  consideration. 

In  the  later  discussion  of  statistical  energy  analysis,  we  can  place  an  analyt¬ 
ical  interpretation  on  0.  In  Eldred’s  approach,  however,  the  emphasis  is  on  the 
empirical  determination  of  0  through  the  analysis  of  missile  vibration  data.  The 
efficiency  of  structural  excitation  depends  on  the  ratio  of  the  acoustic  to  struc¬ 
tural  wavelength.  If  a  structure-  has  all  dimensions  scaled  by  a  constant  factor  K, 
then,  assuming  similarity  in  the  excitation,  the  motion  of  the  original  structure 
at  frequency  /  is  geometrically  preserved  in  the  model  at  frequency  f/K.  Since 
the  acoustic  field  will  preserve  its  similarity  in  scale  K  at  frequency  f/K,  we  can 
expect  that  the  factor  0  will  be  unchanged.  Since  the  radius  of  the  vehicle  a  can 
be  taken  as  a  measure  of  its  size,  Eldied  suggested  that  0  be  considered  a  func¬ 
tion  of  ka  =  2t ra/Xa,  where  \a  is  the  acoustic  wavelength. 
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Since  fSp(f )  is  proportional  to  the  mean  square  pressure  in  constant  per¬ 
centage  bandwidths,  Eq.  (3.6)  would  suggest  that  the  ratio  fSp(f)j{a 2>  should 
depend  only  on  ka.  In  Fig.  14,  we  show  the  values  of  0  as  determined  by  Eldred 
from  rocket  vehicle  vibration  data.  The  data  show  a  spreaj  of  approximately  12 
dB  and  a  general  drop-off  of  5  dB/decade  with  ka.  The  data  are  generally  distrib¬ 
uted  about  (3  ~  1 ,  which  is  the  simple  resonator  value. 

STANOARO  DEVIATION  =  5.88  dB 


Fig.  14.  Experimentally  derived  values  of  factor  P  for  various  rocket  vehicles  [42], 


It  would  be  interesting  to  compare  the  data  of  Eldred  and  Franken  by  re¬ 
ducing  both  to  the  same  form.  Unfortunately,  this  is  made  more  difficult  by 
the  presumption  of  a  single  resonance  in  Eldred’s  approach.  The  original 
measurement  bandwidth  has  been  factored  out  of  Eldred’s  data,  so  that  the 
original  band  pressure  level  corresponding  to  the  measured  value  of  (a2)  cannot 
be  determined.  One  comparison  can  be  made,  however.  Most  of  Franken’s  data 
shows  a  maximum  in  acceptance  or  response  at  2 ra/Ag  —  1,  where  hf  is  the 
longitudinal  wavelength  in  the  structural  material.  This  corresponds  to  ka  = 
\tl\a  ~  15  for  aluminum  or  steel  structures  in  air.  Eldred’s  data  do  not  par¬ 
ticularly  indicate  any  maximum  or  break  in  slope  at  this  value  of  ka.  We  shall 
want  to  examine  implications  of  the  methods  of  both  Franken  and  Eldred  further 
in  the  subsequent  discussions  of  rational  procedures  for  response  estimation. 

3.3  Rational  Approaches  to  Response  Estimation 

What  we  have  termed  rational  approaches  to  response  estimation  are  methods 
of  computation  based  on  the  dynamic  equations  of  motion  of  the  structure. 


ESTIMATION  OF  VEHICLE  RESPONSE 


29 


This  does  not  mean,  of  course,  that  empirical  methods  of  the  previous  section 
are  to  be  regarded  as  irrational.  We  also  note  that  considerable  empiricism  may 
be  involved  in  rational  procedures  in  the  choice  of  structural  models,  selection  of 
parameters,  and  assumptions  regarding  the  loading  by  the  environment. 

The  general  problem  is,  of  course,  the  excitation  of  continuous  structures  by 
loads  that  are  random  in  time  and  are  either  spatially  coherent  (oscillating  shock) 
or  incoherent  (TBL,  BPF,  or  sound).  The  most  direct  approach  is  to  assume  a 
very  ideal  structure,  for  which  the  modes  of  vibration  are  well  known,  and  let  it 
be  excited  by  a  noise  field  of  known  coherence  properties.  This  was  the  approach 
taken  by  Lyon  for  strings  [44, 45] ,  by  Powell  for  membranes  [46] ,  and  by  Dyer 
for  plates  [47] .  In  addition,  Ribner  has  studied  infinite  plate  response  [48]  and 
Cottis  derived  results  for  the  modes  of  a  cylinder  [49] . 

As  helpful  as  much  of  this  work  was,  it  was  recognized  that  actual  structures 
were  more  complex  in  their  behavior  than  the  simple  models  that  had  been 
treated.  Several  attempts  to  treat  more  complex  systems  have  been  made.  They 
include  studies  of  multiple  panel  bays  and  panels  and  beams  with  finite  termin¬ 
ating  impedances  [50-53]. 

Finally,  there  have  been  analytical  studies  that  have  attempted  to  derive 
theoretical  results  in  a  form  directly  comparable  to  the  empirical  methods.  In 
this  approach,  known  generally  as  statistical  energy  analysis,  the  desire  is  to 
predict  average  response  in  frequency  bands.  These  same  techniques  have  also 
been  applied  to  the  transmission  of  vibrational  and  acoustic  energy  to  regions 
within  a  vehicle  [54] . 

Vibration  Estimation  in  Simple  Structures 

The  earliest  studies  of  the  random  excitation  of  continuous  structures  by 
noise  fields  were  concerned  with  Brownian  motion  of  strings,  beams,  etc.,  when 
immersed  in  a  fluid  at  temperature  T.  The  forcing  function  from  the  fluid  is 
generally  taken  to  be  purely  random  in  space  and  time,  with  a  correlation 

ipix,  t )  p(x\  t'))  =  D  5  (x  -  x')  5  (r  -  .'').  (3.7) 

The  parameter  D  is  fixed  by  requiring  that  the  energy  of  vibration  reach  a  value 
of  1/2  kT  for  each  structural  degree  of  freedom  [55] . 

There  are  some  interesting  engineering  problems  for  which  the  loading  corre¬ 
lation  is  similar  to  that  in  Eq.  (3.7).  A  popular  model  for  some  processes  is  rain 
on  the  roof,  short  period  impulses  occurring  randomly  in  time  and  space.  The 
turbulent  boundary  layer,  when  convected  slowly,  may  be  modeled  by  such  a 
process. 

To  demonstrate  the  calculation  procedure,  let  us  imagine  that  the  equation  of 
motion  of  the  responding  system  may  be  written  in  the  form 


(3.8) 
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where  L  is  a  linear  operator  on  the  displacement  w  and  represents  the  elastic 
restoring  forces  in  the  system.  The  damping  is  introduced  through  the  viscous 
coefficient  (3,  and  the  appropriate  mass  density  is  p.  While  this  is  not  the  most 
general  system  one  can  postulate,  it  is  sufficient  for  a  discussion  of  most  methods 
developed  to  date. 

If  the  boundary  conditions  are  specified,  then  the  eigenfunctions  (natural 
modes)  of  the  operator/,  are  determined  from 


L^m  ~  ufa'l'M’  C3-9) 

with  the  normalization  condition 

=  >  (3-10) 

where  8mn  is  the  Kronecker  delta,  and  the  average  is  spatial  average  over  the 
surface  of  the  structure. 

If  the  displacement  u>  and  the  pressure  p  are  expanded  in  a  series  of  the  modal 
shape  functions, 


P(x,t)  =  E  hiWW) 

u<A',  0  =  E  hlW  »W')J 


then  Eq.  (3.8)  becomes 


dt2 


+  p  ~df +  =  plpM’ 


(3.11) 


(3.12) 


where 


PAI(0  =  <p(x,0tM(x))s  1  (313) 

VO  =  M*.')  'l'Af(x)>s  J 

There  are  several  ways  to  obtain  a  solution  of  Eq.  (3.12),  but  a  popular  one  is 
to  use  the  impulse  function  (56) : 

hM(t)  =  —7-  sin  uMt,  (3.14) 


where 


-(«&-* W". 
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If  the  damping  is  small,  ioM  ~  u>’M. 

In  terms  of  the  impulse  response, 

HfoO)  =  P"1  [  dt"hM{t-t")PM(t").  (3.15) 

*^-oo 

The  modal  amplitude  I%(t)  is  of  little  value  in  itself,  since  P^(t')  and  hence 
WM(t)  are  random  functions  of  time  and  are  only  specified  by  their  statistics. 
The  crosscorrelation  of  the  modal  amplitudes  can  be  written 

<wMwN(t'))=p-2  f  dt-  f  dt'"hM(t -  t-)hN(t -  t'ypM(t"wN(n)cns, 

J-OO  J-OO 

(3.16) 

where  the  notation  < .  .  .  .)ens  means  ensemble  average  of  as  distinct  from  the 
spatial  averages  defined  previously.  By  using  Eq.  (3.13),  the  source  term  in 
Eq.  (3.16)  is 

(PM(t")PN(t"'))  =  «^/W^v(A')<p(x,t")p(A,,/"))e/„V>„  (3. 1 7) 

where 

<p(x,  t")p(x',  t"'))cns  =  (p2)  <pp(x  -  x',  t"  -  t'")  (3.18) 

is  the  space-time  pressure  correlation  field.  We  have  discussed  the  form  of  this 
function  for  TBL’s  and  regions  of  separated  flow  in  the  preceding  chapter. 

Using  the  preceding  formulation  and  the  model  of  the  pressure  field  correla¬ 
tion  of  turbulent  flow, 


<{,p(x,t)  =  A,  8(x-Uct)e^°,  (3.19) 

Lyon  [44]  calculated  the  response  of  a  ribbon  to  a  turbulent  pressure  field 
convected  along  its  length.  In  a  similar  set  of  calculations,  Dyer  [47]  computed 
the  response  of  a  rectangular  supported  plate  to  a  simiiar  convected  and  decaying 
correlation  field. 

These  calculations  showed  several  interesting  features  of  the  response  of  struc¬ 
tural  modes  to  a  convected  TBL:  (a)  if  the  convection  speed  is  very  low,  the 
excitation  is  similar  to  rain  on  the  roof;  and  (b)  there  is  augmented  response  for 
modes  that  arc  resonant  and  have  a  trace  wave  speed  in  the  direction  of  con¬ 
vection  that  equals  the  convection  speed.  For  a  given  frequency  this  can  occur 
only  when  the  convection  speed  exceeds  the  speed  of  flexural  waves  on  the 
structure.  The  modal  amplitude  as  a  function  of  convection  speed  is  shown  in 
Fig.  15.  As  we  shall  see,  these  are  the  important  forms  of  response  behavior  for 
structures  that  arc  more  complex  than  simple  strings  or  plates. 
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Fig.  15.  Mean  square  modal  displacement  of  string  as  function  of  convec¬ 
tion  velocity  [45]. 


Let  us  now  assume  that  the  excitation  is  confined  to  a  single  frequency  w. 
We  use  the  complex  convection  Then 


W  -  1  ,  (3.20, 

I'M')  *  IV"'"'  J 

where  hi  and  Ifyj  are  complex  amplitude  factors.  Placing  these  in  Eq.  (3.12), 


Km  = 


~h,lp 


w2(l  +  ip) -u}i' 


(3-21) 


The  mean  square  displacement  in  time  is 
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V-1 


(w2(x,t))t=jRe  iiM{x)i>N(x) 
M,N 


pmpnIp 2 


{co2(l  +  if 3)  -  o4}{w2(l  +  ip)  ■ 


“n}‘ 

(3.22) 


We  can  generalize  the  result  of  Hq.  (3.16)  by  considering  the  Fourier  trans¬ 
form  of  the  pressure: 


p(x,co)  =  lim 


eiojt  dt. 


(3.23) 


If  the  correlation  in  Eq.  (3.18)  depends  only  on  the  time  difference,  then  it  may 
also  be  written 


<p(x',t)p(x",t'))t=j^j  du>P{x,x"\uyii0t,  (3.24) 

where 

P(x',x"-,  co)=  lim  (2T)-1  (p(x,io)p*{x",u))ens  (3.25) 

7— >°° 


is  called  the  cross-power  spectral  density.  If  we  place  Eq.  (3.28)  into  Eq.  (3.13), 
and  substitute  into  Eq.  (3.22),  we  get 


1  00 
(w2(x,t))=—fRe  Y]  ^A/(*)fyv(*)  dco 
p  M,N 


(J\x,x  \(x>)\1>m(x)\Ijn(x  )> 


{p2)  V1  ,  /  \  i  ,  ,  r  /MN(u)do) 

=  L  — i  > 

A/,iV  U  .If  y Af  I 


Iv; 


(3.26) 


where 


and 


yM  =  w  (1  +  iP)  ~  u>jl 


ihd03)  s  <p2  >’ 1  (l\x', ,v";  w) ^/(.t') 'pN(x"))sy  (3-27) 


is  called  the  /owi/  acceptance  of  the  modes  M,N.  This  function  was  originally 
introduced  by  Powell  [46] . 
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The  result  (Eq.  3.26)  has  been  widely  applied  by  design  engineers  because  it  is 
assumed  that  the  mode  shapes,  the  modal  admittances,  and  the  joint  acceptances 
are  readily  computed  or  measured.  However,  particularly  at  higher  frequencies, 
the  excitation  of  response  may  be  dominated  by  rather  subtle  alterations  in  the 
mode  shapes  or  the  spatial  coherence  of  the  excitation. 

The  joint  acceptance  method  has  been  applied  by  Powell  to  strings  [57]  and 
by  his  associates  to  membranes,  panels,  and  other  structures.  In  some  of  his 
papers,  Powell  also  suggested  how  his  methods  might  be  extended  in  an  approxi¬ 
mate  fashion  to  apply  to  more  complex  vibrational  systems  at  higher  frequencies 
[58] .  We  shall  return  to  this  topic  below. 

Applications  to  Complex  Structures 

When  one  tries  to  apply  the  formalism  of  Eqs.  (3.16)  or  (3.26)  to  a  more  com¬ 
plex  structure  such  as  a  flight  vehicle,  then  the  situation  naturally  becomes  much 
more  complex.  The  definition  of  the  crosscorrelation  (or  cross-power  spectrum) 
over  an  extended  surface  is  tedious,  and  the  mode  shapes  for  realistic  structural 
segments  may  not  be  known.  The  determination  of  exact  resonance  frequencies 
and  values  of  modal  loss  factors  a/e  also  difficult,  although  their  average  proper¬ 
ties  can  be  estimated. 

Most  effort  in  ihe  direct  extension  of  the  basic  formulation  has  been  centered 
on  analyzing  the  modal  patterns  of  the  structure.  The  multiply  supported  panel 
has  been  analyzed  by  P.  W.  Smith,  Jr.  [52],  by  Y.  K.  Lin  at  University  of  Illinois 
[50, 51] ,  and  by  Colin  Mercer  at  Southampton  [53] .  The  motivation,  of  course, 
is  the  similarity  between  the  multiply  supported  panel  and  the  multiple  bay 
panels  on  aircraft  and  other  structures.  It  is  found  that  there  is  a  theoretical 
difference  between  the  behavior  of  N  disconnected  panels  and  N  connected 
panels,  but  the  analyses  and  experiments  have  not  been  able  to  demonstrate  that 
the  rather  idealized  connected  panels  arc  more  representative  of  flight  structures 
titan  are  the  simpler  unconnected  panels. 

In  many  ways,  the  large,  modern,  high-speed  digital  computer  seems  to  offer  a 
panacea  for  the  problem  of  structural  complexity.  The  computer  can  be  quite 
helpful  in  calculating  the  lower  order  modes  of  rather  complicated  structures 
consisting  of  combinations  of  stringers,  ring  frames,  reinforced  cutouts,  etc.  The 
mode  shapes  and  resonance  frequencies  arc  obtained  directly  in  the  output.  The 
required  input  is  generally  the  damping,  the  exciting  forces,  and  considerable 
structural  detail.  The  amount  of  time  required  io  program,  perform,  and  process 
the  calculations  depends  in  a  fairly  direct  way  on  the  number  of  calculations  in¬ 
volved.  Th'S  can  easily  be  estimated  by  considering  how  fine-grained  the  com¬ 
putation  must  be  in  frequency  and  space  for  adequate  accuracy. 

Consider  a  1/8-in.  thick  aluminum  shroud,  6  ft  in  diameter  and  15  ft  long. 
Suppose  that  we  arc  interested  in  its  vibratory  response  from  100  to  2500  Hz 
(a  typical  range  of  interest).  The  finest  linear  spatial  scale  is  about  one-quarter 
of  a  bending  wavelength  at  the  highest  frequency  of  interest.  In  our  example. 
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the  wavelength  is  4  in.  The  number  of  mesh-points  is  approximately  4  X  104. 
If  only  a  single  component  of  vibration  (i.e.,  transverse)  is  important  then  this 
represents  the  number  of  simultaneous  equations  to  be  solved  at  each  frequency. 

To  obtain  a  good  statistic  in  the  frequency  domain,  we  require  at  least  one 
calculation  every  modal  bandwidth.  For  the  example,  if  we  choose  a  loss  factor 
of  0.01,  then  at  1000  Hz,  the  bandwidth  is  10  Hz  and  we  would  require  (at 
least)  240  repetitions  of  the  calculations  for  reasonable  spectral  definition. 

The  reader  may  object  that  there  are  interesting  problems  to  be  solved  where 
the  number  of  required  mesh  points  and  frequencies  is  not  so  great.  This  is  true. 
In  the  case  of  sound  and  vibration  transmission,  however,  we  must  include  the 
descriptions  of  acoustic  media  and  connecting  structures.  In  addition,  we  fre¬ 
quently  want  to  know  response  sensitivity  to  changes  in  parameters  such  as 
damping,  skin  thickness,  and  frame  spacing.  Even  in  modest  systems,  the  number 
of  degrees  of  freedom  increases  rapidly  as  we  include  the  frequencies  of  direct 
engineering  interest. 

Because  of  the  above  considerations  and  others  such  as  cost,  time,  and  output 
format,  digital  computer  calculations  using  finite  element  models  have  been  of 
quite  marginal  value  in  analyzing  sound  and  vibration  response  of  structures. 
There  is  considerable  literature  on  the  subject,  and  we  feel  that  including  it 
hero  would  dilute  the  material  that  we  regard  as  more  significant  to  our  present 
purposes. 

Beginning  in  the  early  1960’s,  a  new  approach  to  the  excitation  of  complex 
structures  by  random  environments  began  to  appear.  These  studies,  known  col¬ 
lectively  as  statistical  energ)’  analysis  (SEA)  have  been  centered  mostly  at  Bolt 
Bcranek  and  Newman,  but  other  groups  have  also  made  contributions  to  the 
literature  in  this  area.  Broadly  speaking,  an  attempt  is  made  to  model  the  struc¬ 
tural  system  by  an  assemblage  of  modes  whose  primary  parameters  (resonance 
frequency,  damping,  mode  shape,  etc.)  are  statistically  determined.  The  primary 
dynamic  variable  is  chosen  to  be  the  energy  of  vibration,  from  which  other 
variables  (acceleration,  strain,  etc.)  may  be  determined.  Much  of  the  inspiration 
for  this  point  of  view  has  been  derived  from  the  methods  of  the  statistical 
analysis  of  room  acoustics. 

Statistical  modeling  and  the  use  of  energy  variables  would  likely  be  valuable 
additions  to  the  too!  kit  of  the  analyst  in  themselves,  but  two  other  principles 
increase  the  power  of  SEA  appreciably,  it  has  been  shown  that  parameter  sta¬ 
tistics  calculated  for  quite  simple  systems  may  be  applied  to  more  complex  sys¬ 
tems.  For  example,  the  average  modal  density  (number  of  modal  resonances 
occurring  in  a  frequency  interval)  may  be  computed  for  the  flexural  resonances 
of  a  simple  rectangular  supported  panel,  but  the  result  of  this  calculation  is  ap¬ 
plicable  to  any  panel  of  the  same  area,  when  one  considers  resonance  frequencies 
wcli  above  the  fundamental. 

Tire  second  principle  states  that  under  random  excitation,  the  flow  of  me¬ 
chanical  power  from  one  system  mode  to  anoti.er  is  proportional  to  the  dif¬ 
ference  in  their  time  average  energies.  This  principle  allows  one  to  infer  the 
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direction  and  magnitude  of  energy  flow  in  a  system,  once  the  input  power  and 
dissipative  processes  are  determined.  Such  calculations  are  necessary  when  con¬ 
sidering  the  transmission  of  sound  and  vibration  in  flight  vehicles. 

It  would  be  very  difficult  to  review  all  the  applications  that  have  been  made 
of  SEA.  Rather  than  do  this,  we  shall  summarize  briefly  many  of  these  calcula¬ 
tions  and  experiments  in  the  following  few  paragraphs.  Subsequent  to  this,  we 
shall  review  in  rather  more  detail  its  application  to  a  spacecraft-shroud  combina¬ 
tion  excited  by  a  sound  field. 

Review  of  SEA  Applications 

The  earliest  studies  that  applied  the  ideas  of  modal  density  and  asymptotic 
behavior  of  mechanical  systems  to  structural  response  must  include  those  of 
Powell  [58]  and  Skudrzyk  [59]  in  1958.  Powell  showed  that  certain  response 
averages  over  modes  were  equivalent  to  the  same  averages  over  infinite  structures. 
In  his  analysis,  Skudrzyk  examined  the  input  impedance  of  multimodal  struc¬ 
tures  and  showed  that  the  modal  density  has  a  central  role  in  the  bridge  between 
asymptotic  modal  behavior  and  behavior  of  the  same  system  infinitely  extended. 

The  response  of  infinite  plates  and  cylinders  to  a  TBL  pressure  field  were 
carried  out  by  Corcos  and  Liepmann  [60]  and  by  Ribner  [48] .  Dyer,  however, 
appears  to  be  the  first  to  group  resonant  modal  response  into  frequency  bands, 
and  develop  response  spectra  for  finite  systems.  He  used  the  calculations  of 
modal  response  to  TBL  pressures,  and  the  modal  density  of  a  flat  plate  to  pre¬ 
dict  the  displacement  spectrum.  The  result  that  Dyer  obtained  is  fully  equivalent 
to  the  infinite  system  results  of  Corcos  and  Liepmann. 

The  panel  response  to  a  TBL  pressure  field  is  simplified  if  one  does  not  con¬ 
sider  the  reaction  of  the  fluid  to  panel  vibration  (except  in  average  ways  rep¬ 
resented  by  mass  loading  and  radiation  damping).  The  response  of  a  panel  to 
sound  is  very  much  related  to  its  ability  to  radiate  sound.  Smith  has  shown  that 
reciprocity  arguments  can  be  applied  to  acoustically  induced  panel  response 
calculations  [61].  Since  the  energy  exchange  between  the  sound  field  and  the 
structure  is  conservative,  and  both  have  energy  storage  capabilities,  the  need  for 
a  principle  of  energy  sharing  between  randomly  excited  systems  is  evident. 

The  simplest  system  for  the  study  of  energy  sharing  between  mechanical  sys¬ 
tems  is  a  pair  of  simple  resonators.  The  power  flow  between  two  randomly  ex¬ 
cited  resonators  with  weak,  linear  coupling  was  studied  by  Lyon  and  Maidanik 
[62] .  The  result  showed  that  the  power  flow  was  proportional  to  the  difference 
in  time  average  energy  of  the  resonators  and  that  it  flows  from  the  resonator  of 
higher  energy  to  the  one  of  lower  energy.  More  recently,  Scharton  [63]  and 
Ungar  [64]  have  shown  that  the  weak  coupling  restriction  can  be  removed,  and 
Newland  [65]  has  extended  the  result  to  include  certain  forms  of  nonlinear 
coupling  elements. 

The  two-resonator  system  is  of  relatively  little  interest  in  itself.  One  is  usually 
concerned  about  sets  of  resonators  (modes)  that  may  exchange  energy,  if  one 
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can  argue  (or  demonstrate)  that  the  energies  of  the  modes  in  a  particular  group 
are  nearly  equal,  then  the  interactions  of  that  modal  group  with  other  mode  sets 
can  be  calculated  using  the  single  resonator  pair  results.  This  approach  is  taken 
by  Lyon  [66] ,  Maidanik  [67] ,  and  Manning  and  Maidanik  [68]  in  studying  the 
interaction  of  an  elastic  panel  with  sound.  The  energy  interaction  between  two 
elastic  structures  has  been  studied  by  Lyon  and  Eichler  [69]  and  by  Scharton 
[63] .  A  recent  survey  of  the  relation  between  multimodal  and  mode-pair  energy 
sharing  has  been  prepared  by  Ungar  [64] . 

Finally,  there  is  a  set  of  problems  concerning  the  transmission  of  energy 
through  a  system  by  one  or  more  paths.  The  system  that  we  shall  discuss  in 
some  detail  below  is  of  this  type.  Studies  of  sound  transmission  through  struc¬ 
tures  using  SEA  have  been  made  by  Lyon  [70],  Eichler  [71]  and  White  [72]. 
A  study  of  vibration  transmission  in  a  3-element  structure  has  been  carried  out  by 
Lyon  and  Scharton  [73] .  A  study  of  sound  and  vibration  transmission  in  a 
spacecraft  system  has  been  completed  recently  by  Manning,  Lyon  and  Scharton 
[7] .  This  is  the  example  that  we  shall  treat  below. 

In  summary,  the  SEA  treats  energy  distributions  in  multi-element  systems  by 
(a)  computing  mechanical  power  injected  into  the  elements  by  the  environment, 
and  (b)  by  calculating  the  redistribution  of  this  energy  throughout  the  system 
due  to  dissipative  processes  and  to  inter-element  coupling.  It  is  still  a  set  of 
procedures  under  development,  but  its  results  to  date  suggest  that  it  has  an  im¬ 
portant  role  to  play  in  prediction  of  response  to  environmental  loads. 

Example  of  Spacecraft  Vibration  Prediction 

The  vibration  of  a  spacecraft  at  launch  is  primarily  due  to  the  acoustic  noise 
exterior  to  the  vehicle.  Let  us  considei  the  spacecraft  shown  in  Fig.  16  as  an 
example.  Since  a  sound  field  stores  mechanical  energy,  we  call  this  element  1. 
The  energy  of  the  sound  field  is  converted  into  vibrational  energy  by  the  ex¬ 
terior  structure  or  shroud  of  the  spacecraft.  This  is  element  2.  The  exterior 
structure  in  turn  radiates  acoustic  energy  into  the  air  space  contained  within  the 
structure.  This  air  space  is  element  3.  Finally,  the  structural  panels  of  the  space¬ 
craft,  element  4,  are  excited  by  this  sound  field  into  vibration. 

Figure  16  shows  an  additional  path  of  energy  transmission.  This  path  includes 
the  ring  frame,  clement  5,  the  mounting  trusses,  element  6,  and  the  spacecraft. 
We  shall  want  to  include  these  elements  as  a  path  for  mechanical  energy  in 
parallel  to  the  acoustic  path.  These  paihs  are  shown  diagramatically  in  Fig.  17. 

Consider  the  acoustic  path  first.  We  imagine  the  exterior  sound  field  to  be 
contained  in  a  large  room  of  volume  V.  If  the  space-time  rms  pressure  in  this 
field  in  a  band  of  width  A  is  px ,  then  the  average  energy  of  an  acoustic  mode  is 


?j]A  pc 2  «jA 


(3.28) 
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The  first  factor  is  the  acoustic  energy  density  and  nx  is  the  number  of  acoustic 
modes  per  unit  frequency  bandwidth.  This  quantity  is  termed  the  modal  density 
and  is  given  by 


n,(0“  4tt2/2T7c3.  (3.29) 

Note  that  the  volume  V cancels  out  of  Eq.  (3.29). 

The  energy  balance  equation  for  the  shroud  must  include  the  power  supplied 
to  it  by  the  sound  field,  the  power  that  it  radiates  to  the  interior  space,  and  the 
power  that  it  dissipates  internally.  Formally 

III, a  =  n2(3  +  nf (3.30) 

or,  introducing  dissipative  and  coupling  loss-factors, 
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(3.31) 


The  left  side  of  this  equation  shows  the  proportionality  of  power  flow  from  1  to 
2  on  the  difference  in  modal  energies  of  sound  field  and  shroud.  The  first  term 
on  the  right  is  the  power  flow  from  the  shroud  to  the  interior  space,  while  the 
final  term  represents  the  dissipation  in  the  shroud.  The  quantities  77 1>2  andi72i3 
are  the  “coupling  loss  factors,”  and  they  must  be  calculated  or  measured  for  the 
system.  The  dissipation  loss  factor  r?2  is  usually  obtained  experimentally. 

If  the  transmission  of  sound  to  element  3  is  temporarily  ignored,  then 
t?2> 3  0  in  Eq.  (3.31)  and  the  ratio  E2/Ei  can  be  obtained  from  Eq.  (3.31). 
The  modal  density  n2  is  derived  from  shell  dynamics  of  cylinders  [74,  75]  and 
r?12  is  obtained  by  using  a  consistency  relationship. 

The  left-hand  side  of  Eq.  (3.31)  represents  the  power  flow  from  element  1  to 
clement  2.  If  we  were  writing  an  energy  balance  for  clement  1,  we  would  include 
this  same  term  in  the  form 


-W2,i>h 


since  it  represents  energy  lost  from  element  1.  But  these  power  flows  are  equal. 
We  therefore  have  the  consistency  relationship 

=  na,ih.i.  (3.32). 

The  importance  of  this  is  that  the  modal  densities  are  relatively  easy  to  calculate, 
and  if  we  know  one  of  the  coupling  loss  factors,  we  do  not  have  to  compute  the 
other.  In  this  particular  instance,  172,1  is  known  since  it  is  very  simply  related  to 
the  radiation  resistance  of  a  cylindrical  shell  [67] : 
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The  modal  density  n2  and  coupling  loss  factor  are  shown  in  Figs.  18  and  19  for 
a  glass  fiber  shroud  195  in.  long,  65  in.  in  diameter  and  0.13  in.  thick.  Note  the 
peak  in  both  the  modal  density  and  the  coupling  loss  factor  at  the  ring  frequency. 
The  energy  E2  is  related  to  the  space-time  mean  square  transverse  acceleration  of 
the  cylinder  by 


E2  =  M2(a\)!o?. 


(3.34) 


Fig.  18.  Modal  density  of  cylinder,  averaged  in  one-third  octave 
bands,  compared  to  fiat  plate  of  equal  area  (7). 


If  in  addition,  we  assume  loss  factors  of  y2  =  10"2  and  3.16  X  10-2  for  the 
shroud,  then  the  predicted  acoustically  induced  response  is  shown  in  Fig.  20. 
It  is  interesting  to  compare  this  calculated  response  with  the  semi-empirical  pre¬ 
diction  of  Frankcn  in  Fig.  13. 
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THIRO  OCTAVE  BAND  CENTER  FREQUENCY  ( Hz > 

Fig.  19.  Shroud-sound  field  coupling  loss  factor  and  dissipative  loss  factor  of  shroud  [7]. 


Fig.  20.  Acoustic  acceptance  of  shroud  for  two  values  of  structural  damping  [7] . 
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The  energy  equations  for  the  interior  acoustic  space  are  similar  to  Eq.  (3.30): 

n2(3  +  n^  =  n3>4  +  nf **  (3.35) 

or,  again  introducing  the  loss  factors, 


f E2  £3  \  _  (E3  £<■  \  v  -3/r\ 

w«2??2,3  — - +  nt  f  -  wn3i?3 ,4  — - +  V3E3.  (3.36) 

\«2  «3  /  \«3  «4  / 

In  this  equation,  an  additional  term  has  been  added  to  account  for  power  flow 
directly  from  the  external  sound  field  to  the  enclosed  space  due  to  the  vibration 
of  shroud  modes  that  resonate  outside  the  frequency  band  A.  In  ordinary  acous¬ 
tics,  this  is  the  mass  law  contribution  to  sound  transmission  [76] .  Since  this 
power  flow  is  again  proportional  to  the  differences  in  modal  energies  of  systems 
1  and  3, 


nNR  _  (Ei  E3\ 
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(3.37) 


where  the  coupling  loss  factor  th>3  is  related  to  the  acoustic  transmission  loss  of 
the  shroud. 

Finally,  the  acoustically  induced  energy  of  the  spacecraft  E 4  is  determined  by 
its  power  balance  relations 


n  =  ndi5S 

ll3,4  1  *4  > 


or,  using  loss  factors  and  energies, 

\«3  «4  / 
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(3.38) 


(3.39) 


since  all  the  energy  supplied  to  the  spacecraft  by  the  sound  field  in  element  3 
must  be  dissipated  (neglecting  the  connection  to  the  truss).  Relating  the  energy 
£4  to  the  surface  acceleration  by 


Ea  =  MA(al)lo}2 


(3.40) 


and  then  by  using  Eqs.  (3.31),  (3.36)  and  (3.39),  we  can  generate  a  ratio  of  ac¬ 
celeration  a4  to  external  sound  pressure  px .  This  ratio  is  shown  logarithmically 
in  Fig.  21  as  acoustic  path  only. 

In  calculating  the  acoustic  transmission,  the  mechanical  path  through  the  ring 
frame  (5)  and  truss  (6)  were  neglected.  The  mechanical  power  supplied  to  the 
spacecraft  by  the  truss  is  derived  from  the  following  power  balance  relations: 
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Fig.  21.  Vibrational  acceptance  of  spacecraft  due  to  acoustic  and  vibrational  paths 
separately  (7] . 


1.  For  element  5,  the  ring  frame. 
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(3.41) 


(3.42) 


2.  For  element  6,  the  mounting  truss, 


n„  -  nM  +  nf“ 


(3.43) 


or 


%E‘  '■  (3-44) 


3.  For  element  4,  the  spacecraft, 


n  =  n'w  =  nrfte 

u6,4  11 5,4  u4 


(3.45) 
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or 


^6,4 »6  (“  ” ' ~  )  +  n?4  =  <*n4EA  •  (3-46) 

Again  in  Eq.  (3.43)  we  have  included  the  possibility  of  energy  transmission  from 
the  ring  frame  to  the  spacecraft  by  nonresonant  motion  of  the  truss.  Formally, 
this  is  given  by 


n 


NR 

5,4 


(3-47) 


The  calculations  of  the  coupling  loss  factors  are  carried  out  by  Manning  et  al. 
[7]  for  the  structural  elements  shown  in  Figs.  22,  23,  and  24.  The  energy  E2  in 
Eq.(3.39)  is  assumed  known  by  application  of  Fig.  15.  From  Eqs.  (3.42),  (3.44) 
and  (3.46),  one  can  calculate  the  ratio  of  spacecraft  acceleration  to  external 
pressure  due  to  the  mechanical  paths.  This  has  been  done  and  is  plotted  in 
Fig.  21  as  mechanical  path  only. 
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Fig.  22.  Channel  beam  model  of  ring  frame,  element  5  [7] . 


In  this  example,  the  result  of  the  calculations  showed  that  the  acoustic  path 
was  dominant,  in  that  it  produced  spacecraft  vibration  levels  of  the  order  of  10 
to  20  dB  greater  than  did  the  mechanical  path.  We  hasten  to  point  out  that  this 
is  a  particular  result  depending  on  the  geometry  and  structural  parameters  of  the 
chosen  configuration.  Other  spacecraft  in  a  different  shroud  may  have  signifi¬ 
cantly  different  relative  sensitivities  to  the  acoustic  and  mechanical  paths. 

In  discussing  the  use  of  SEA  in  this  example,  we  have  only  been  able  to  indi¬ 
cate  the  nature  and  order  of  the  calculations.  Much  additional  detail  can  be  found 
in  Ref.  7.  It  should  be  clear  from  the  above,  however,  that  while  SEA  offers 
substantia!  simplifications  in  vibration  prediction  over  “classical”  methods,  it  still 
may  involve  relatively  lengthy  computations.  The  length  and  difficulty  of  the 


Fig.  24.  Diagram  of  construction  of  spacecraft  panels,  element  4  [77] 
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calculations  is  usually  dependent  on  how  much  cleverness  and  intuition  one  can 
bring  to  the  analysis.  As  additional  experience  is  gained  in  these  studies,  sim¬ 
plifications  in  the  analyses  will  follow.  It  is  unlikely  that  they  will  ever  achieve 
the  simplicity  of  the  semi-empirical  procedures,  however,  but  the  range  of  prob¬ 
lems  that  SEA  is  able  to  handle  is  potentially  much  greater. 

3.4  Conclusions 

The  prediction  of  vehicle  response  to  fluctuating  loads  is  still  in  a  primitive 
state  as  far  as  the  designer  or  test  engineer  is  concerned.  Generally,  the  purely 
empirical  methods  are  currently  used  for  response  prediction,  although  some  of 
the  semi-empirical  procedures  may  be  used  on  existing  data  to  scale  for  changes 
in  vehicle  diameter  and  surface  mass  density.  The  use  of  rational  procedures, 
such  as  SEA,  is  still  in  the  research  stages. 

An  important  method  of  determining  vibratory  response  to  loads  is  the  use  of 
experiments  This  method  simulates  flight  loads,  and  response  levels  are  inferred 
by  direct  measurement.  Since  such  studies  are  very  close  to  the  testing  proce¬ 
dures  to  be  discussed  in  the  following  chapter,  we  shall  defer  discussing  them. 


CHAPTER  4 

TESTING  FOR  ENVIRONMENTAL  EFFECTS 


The  role  of  the  laboratory  in  environmental  studies  is  very  important.  Ulti¬ 
mately,  the  best  judgment  of  the  adequacy  of  a  design  to  withstand  the  environ¬ 
ment  must  rely  on  laboratory  acceptance  or  proof  tests  of  the  system  and/or  its 
components.  Well-chosen  simulation  of  the  service  environment  will  optimize 
the  chances  for  a  successful  system,  while  avoiding  the  excessive  time  and  costs 
of  overdesign. 

The  environmental  laboratory  has  a  second  role,  however,  that  may  become 
obsc’ired  in  the  pressures  of  system  scheduling.  The  establishment  of  important 
structural  parameters,  or  the  determination  of  dominant  modes  of  behavior  as  a 
check  or  substitute  for  theoretical  derivations,  is  an  important  function  for  the 
laboratory.  In  this  chapter,  we  discuss  both  of  these  functions. 

4.1  Purposes  of  Environmental  Testing 

The  traditional  proof  testing  activity  of  the  environmental  laboratory  is 
basically  one  of  environmental  simulation.  The  philosophy  is:  the  equipment 
is  to  be  subjected  to  loads  that  represent  the  most  reasonably  severe  acoustic 
and/or  vibration  environment  likely  to  be  encountered  during  its  mission(s).  The 
duration  of  the  test  is  set  by  an  analysis  of  the  Hm»  hict^r*  of  the  load  intensity, 
either  derived  from  flight  tests  or  from  theoretical  analysis.  In  some  tests,  such 
as  those  for  fatigue  of  structures,  the  time  scale  of  the  test  may  be  reduced  by 
increasing  the  intensity  of  the  excitation  [78] . 

Experimental  evaluation  of  response  parameters  becomes  even  more  important 
when  statistical  energy  analyses  (SEA)  is  used  as  a  framework  for  understanding 
structural  response.  As  an  example,  the  relative  role  of  resonant  and  nonresonant 
vibration  in  overall  response  or  transmission  of  energy  can  be  assessed  by 
changing  the  structural  damping.  Also,  the  importance  of  a  vibration  transmis¬ 
sion  path  vis  a  vis  an  acoustic  path  may  be  determined  by  significantly  changing 
the  characteristics  of  either  path;  c.g.,  by  breaking  the  structural  connection  or 
by  adding  an  absorptive  treatment  to  the  acoustic  spaces. 

The  parameter  evaluation  process  just  described  is  roughly  equivalent  to  what 
has  been  previously  described  as  experimental  analysis  [79] ,  which  is  the  use  of 
experiments  to  ask  fundamental  questions  about  system  behavior.  Experimental 
analysis  is  a  complement  to  theoretical  analysis,  in  that  it  may  suggest  system 
models  that  are  more  amenable  to  theoretical  study.  This  process  is  to  be 
contrasted  to  the  proof  testing  described  earlier  (which  asks  essentially  one 
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question  about  the  equipment-did  it  break?)  and  the  research  tests  described 
below. 

By  research  tests  are  meant  experiments  for  determining  the  basic  acoustic 
and  mechanical  information  on  a  system.  Such  tests  include  the  measurements 
of  resonance  frequencies  and  mode  shapes  [80] ,  the  structural  damping  values 
[81] ,  and  coupling  loss  factors.  The  results  of  these  tests  provide  structural  in¬ 
formation  for  input  to  theoretical  or  experimental  analyses.  We  shall  not  have 
much  to  say  about  research  tests  in  this  chapter,  but  guidelines  for  such  testing 
are  given  in  the  references  [7,  Appendix] . 

The  parameter  evaluation  function  of  the  environmental  laboratory  is 
beginning  to  receive  more  attention  as  semi-empirical  and  rational  methods  of 
response  estimation  become  more  widely  used.  The  methods  of  Franken  [41] 
and  Eldred  [42]  suggest  that  vehicle  dimensions  and  weights  will  normalize  the 
response  to  pressure  ratios.  Laboratory  tests  can  readily  verify  whether  these 
schemes  apply  to  the  vehicle  at  hand. 

4.2  How  Specifications  Are  Set 

In  chapter  3  we  discussed  how  a  particular  vibration  or  acoustic  spectrum 
measured  in  flight  might  be  converted  to  a  vibration  specification.  lr-  this  section, 
we  broach  a  rather  more  complex  question:  how  is  (or  should)  an  environmental 
specification  be  developed  based  on  mission  analysis,  failure  experience,  and  the 
manner  of  simulation.  Let  us  review  how  each  of  these  has  been  used  for 
developing  specifications. 

The  largest  difference  between  missions  is  the  one-shot  life  of  space  boosters, 
and  the  multi-flight  lifetime  of  aircraft.  The  collection  of  time  history  data  for 
each  vehicle  type  is  quite  different.  Rather  complete  mission  profile  data  are 
collected  for  each  spacecraft-booster  combination,  at  least  for  frequency  bands 
included  in  the  telemetry  capabilities  [3J .  Surprisingly,  long  time  histories  are 
not  as  available  for  aircraft,  although  steps  have  recently  been  taken  to  gather 
such  data  for  critical  sections  of  aircraft  [82] . 

The  type  of  mission  profile  has  strong  implications  for  the  kind  of  failure 
mode  that  is  likely  to  be  most  significant.  The  multi-flight  mission,  involving 
many  hours  of  flight  with  intensive  fluctuating  pressure  fields  can  cause  metal 
fatigue  to  occur  in  many  important  structural  segment!  The  service  environ¬ 
ment  may  endure  so  long  that  it  is  impractical  to  test  for  failure  at  expected 
levels.  In  these  cases,  the  levels  may  be  increased  by  a  series  of  factors  and  the 
fatigue  life  recorded  for  each  set  of  levels.  The  lifetime  under  service  loads  is 
then  obtained  by  extrapolation.  The  shape  of  the  excitation  vs  life  curves  is 
usually  such  that  this  will  produce  a  conservative  estimate  of  the  lifetime  at 
lower  excitation  levels  [83] . 

Malfunction  of  single  flight  spacecraft  has  usually  been  identified  with  failure 
of  on-board  equipment.  The  environmental  problem  is  more  closely  related  to 
an  exceedance  of  tolerable  limits  raihei  than  a  fatigue  process.  Accelerated, 
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augmented  amplitude  testing  does  not  make  sense  for  this  kind  of  failure  mode. 
Tests  of  spacecraft  structure  and  components  are  carried  out  at  levels  and  dura¬ 
tions  that  simulate  the  actual  mission  profile  as  closely  as  possible,  a  practicable 
process  for  the  shorter  duration  environment  of  spacecraft. 

The  goal  of  the  proof  test  is  primarily  a  demonstration  of  an  adequate  design. 
Exact  simulation  of  the  loading  (with  a  safety  factor)  is  one  way  to  achieve  this. 
Exact  simulation  is,  however,  almost  always  beyond  the  state  of  the  art  and  is 
always  expensive.  It  is  necessary,  therefore,  to  identify  the  most  critical  part  of 
the  environment  and  the  attendant  loads.  For  example,  in  some  spacecraft,  the 
greatest  vibration  of  electronic  components  may  occur  during  lift-off,  and  in 
other  spacecraft,  this  may  occur  during  transonic  flight.  In  the  former,  the 
vibration  is  acoustically  induced,  while  in  the  latter,  the  excitation  may  be  a 
combination  of  separated  flow  and  oscillating  shock. 

To  continue  with  the  example,  the  component  that  fails  probably  does  so 
because  of  mechanical  strain  or  acceleration  in  its  mounting  and  not  because  of 
the  acoustic  pressure  surrounding  it.  If  it  sits  on  a  circuit  board,  the  story  is 
different-the  circuit  board  may  be  quite  sensitive  to  the  sound  field,  and  the 
acoistic  pressures  surrounding  the  circuit  board  may  be  quite  dependent  on 
whether  the  external  pressure  environment  is  due  to  acoustic  or  separated  flow 
[7,  section  3].  Thus,  simulation  is  seen  to  be  a  sophisticated  concept.  It  might 
be  adequate  to  simulate  the  vibration  of  the  component  or  the  loads  on  the 
spacecraft,  but  a  simple  vibration  or  acoustic  simulation  at  intermediate  points 
in  the  energy  transmission  chain  may  not  be  adequate.  At  this  point,  experi¬ 
mental  analysis  may  come  into  play,  since  our  evaluation  studies  may  have  shown 
that  the  circuit  board  was  much  more  sensitive  to  the  vibration  of  its  mounting 
frame  than  to  expected  sound  levels.  We  would  then  conclude  that  a  vibration 
simulation  for  the  board  mounting  was  acceptable. 

In  the  preceding  paragraphs,  we  have  assumed  that  the  environment  can  be 
known,  and  simulated,  if  necessary,  at  any  point  in  the  energy  transmission  chain. 
As  we  discussed  in  chapter  2,  a  service  noise  environment  is  so  complex  that  one 
is  forced  to  model  it  with  an  environment  achievable  in  the  laboratory.  This 
environment  may  be  a  reverberant  sound  field,  or  three-axis  vibration  excitation 
at  a  mounting  point  for  a  piece  of  equipment. 

A  separate  but  related  problem  to  that  of  detail  is  the  role  of  the  dynamic 
interaction  between  systems.  Thus,  one  may  wonder  if  the  vibration  of  the  sur¬ 
face  of  a  vehicle  will  produce  a  break  reaction  on  the  fluid  to  change  acoustic  or 
aerodynamic  pressures.  Theoretically,  if  one  could  produce  the  service  pressures 
at  the  vehicle  surface  under  conditions  of  service  vibration,  then  the  simulation 
would  be  accomplished.  The  same  comments  apply  to  a  vibrational  problem  in 
which  a  mounting  truss  or  bulkhead  is  mounted  to  a  ring  frame. 

The  difficulty  in  all  this  is  reproduction  of  the  >ervice  levels.  The  dynamic 
interaction  between  the  two  systems  is  embedded  in  the  statistical  covariances  of 
the  dynamic  variables  at  the  system  boundaries.  To  reproduce  the  interaction 
effects  when  imposing  the  environment  at  an  interface,  it  would  be  necessary  to 
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reproduce  these  covariances,  which  fluctuate  rapidly  as  resonances  of  the  test 
structure  are  encountered.  It  is  usually  not  possible  (or  desirable)  to  attempt 
such  a  detailed  simulation  even  if  the  interface  data  were  known  to  such  detail. 
For  this  reason,  another  approach  to  handling  interaction  effects  is  to  simulate 
the  dynamics  at  the  interface  and  the  environmental  levels  at  a  location  away 
from  the  interface. 

A  simple  example  of  this  process  occurs  in  the  acoustic  testing  of  component 
packages  by  diffuse  sound  fields.  At  very  low  frequencies,  the  acoustic  wave¬ 
lengths  are  large,  and  the  component  is  merely  squeezed  by  the  sound  field.  At 
very  high  frequencies,  the  wavelengths  are  small  compared  to  the  dimensions  of 
the  package,  and  one  has  pressure  correlations  on  the  package  that  are  quite 
similar  to  those  in  the  reverberant  field.  In  the  frequency  regime  where  wave¬ 
lengths  are  of  the  order  of  package  dimensions,  however,  acoustic  diffraction 
effects  can  result  in  highly  complicated  spatial  sound  pressure  correlations  that 
are  quite  difficult  to  predict  or  to  specify.  The  usual  method  for  specifying  the 
acoustic  environment  in  such  cases  is  to  require  that  the  sound  field  be  rever¬ 
berant  and  that  its  spectrum  be  known  at  a  location  away  from  the  test  object. 
The  interaction  effects  are  thus  taken  care  of  through  the  dynamics  of  the  sound 
field-component  interaction  during  the  test. 

Although  not  done  as  commonly  as  in  acoustics,  it  is  possible  to  approach  the 
vibra'ion  interaction  problem  in  the  same  way.  We  shall  have  more  to  say  about 
this  in  section  4.3. 

The  definition  of  a  specification  is,  therefore,  a  systems  problem,  involving 
the  use  of  flight  data  analysis  of  missions,  limitations  and  cost  of  facilities,  pos¬ 
sible  failure  modes,  and  desired  reliability.  Rarely  are  all  these  things  taken  into 
account  in  the  environmental  specification  given  to  the  equipment  supplier.  By 
the  time  the  supplier  has  proof  tested  and  delivered  the  item,  however,  most  of 
them  have  had  the  most  agonizing  kind  of  consideration! 

4.3  Direct  Environmental  Simulation 

By  direct  environmental  simulation,  1  mean  the  replacement  of  an  acoustic 
environment  with  a  sound  field,  or  of  a  vibration  environment  by  a  mechanical 
shaker.  Simulations  that  replace  a  TBL  with  a  sound  field,  or  a  sound  field  with 
a  shaker  will  be  discussed  in  the  following  section. 

Acoustic  Simulation 

The  overall  sound  pressure  level  on  flight  structures  may  reach  as  much  as  170 
dB  re  0.0002  jubar.  In  a  progressive  wave,  this  amounts  to  nearly  I04  w?tts/sq  ft. 
For  a  structural  segment  of  reasonable  dimensions,  this  means  that  acoustic 
power  of  the  order  of  104  to  106  watts  is  required  to  provide  the  necessary 
levels.  It  is,  of  course,  very  difficult  and  very  expensive  to  generate  acoustic 
power  of  this  order  of  magnitude. 
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There  are  two  large  acoustic  chambers  in  the  country  designed  to  produce 
sound  levels  of  this  magnitude  on  large  vehicle  segments.  One  of  these  is  the 
Sonic  Fatigue  Facility  [84]  at  Wright  Field  in  Dayton,  Ohio,  shown  in  Fig.  25. 
The  other  is  the  Launch  Environment  Acoustic  Facility  [85]  at  the  Manned 
Spacecraft  Center  in  Houston,  Texas.  The  first  of  these  uses  amplitude  and 
frequency  modulated  sirens  as  sources,  while  the  NASA  facility  uses  air  mod¬ 
ulated  horns.  Both  facilities  are  designed  for  traveling  wave  and  reverberant 
field  modes  of  operation. 


Wright-Pattcrson  Air  Force  Base  photograph. 
Fig.  25.  Sketch  of  Wright  Field  Sonic  Fatigue  Facility. 


The  production  of  a  certain  frequency  spectrum  of  pressure  over  a  range  of 
positions  can  usually  be  achieved  satisfactorily.  Specification  writers  generally 
acknowledge  that  the  complications  of  directivity  of  the  sound  cannot  be 
achieved,  so  the  requirements  are  written  as  reverberant  or  traveling  wave.  The 
spatial  distribution  of  sound  levels  over  the  surface  of  the  structure  is  also  dif¬ 
ficult  to  simulate  in  detail,  but  a  fair  latitude  is  usually  allowed  (in  part  because 
the  environmental  information  is  quite  limited  in  this  regard). 
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The  acoustic  space  may  have  a  rather  subtle  effect  on  panel  dynamics.  If  the 
sound  field  is  concentrated  about  the  structure  by  close  fitting  duct  work,  then 
the  damping  and  resonance  frequencies  of  the  structure  may  be  perturbed.  The 
increased  coupling  may  damp  the  panel  so  strongly  that  the  fundamental  panel 
resonance  does  not  appear  [86]. 

In  special  circumstances,  it  has  been  desirable  to  test  a  vehicle  outdoors.  There 
are  several  reasons  for  this:  (a)  a  semi-anechoic  environment  may  be  desired, 
(b)  the  size  of  the  assembled  vehicle  may  preclude  chamber  testing,  or  (c)  the 
captive  operation  of  other  engines  may  provide  a  convenient  source  of  high 
intensity  sound  for  structural  or  equipment  testing.  Each  advantage,  of  course, 
implies  certain  disadvantages:  (a)  the  directional  properties  of  the  acoustic  field 
are  likely  not  representative  of  service  environment,  (b)  the  provision  of  handling 
equipment  may  be  quite  expensive,  and  (c)  engine  testing  schedules  may  not  be 
convenient  for  the  structural  tests.  Nevertheless,  such  tests  have  proven  quite 
useful  in  the  past  and  offer  an  attractive  and  potentially  inexpensive  alternative 
to  large  facility  construction. 

Finally,  we  note  a  severe  and,  at  present,  fundamental  limitation  of  acoustic 
sources-nonlinearity.  As  we  noted  above,  single  intense  sound  sources  may  be 
expected  to  generate  power  of  the  order  of  103  to  10s  watts  which  is  carried 
by  the  sound  wave.  In  a  1-sq  ft  duct,  these  correspond  to  sound  pressure  levels 
Lp  from  160  to  180  dB  re  0.0002  pbar.  The  progression  of  a  pure  tone  (single 
frequency)  wave  at  these  levels  will  have  the  form  shown  in  Fig.  26.  In  this 
progression,  the  positive  peak  of  the  wave  advances  until  it  becomes  triangular 
after  a  distance  [87] 


xc  =  X/8M,  (4.1) 

where  M  =  u/c  is  the  Mach  number  o'  'he  acoustic  particle  velocity  and  X  is  the 
acoustic  wavelength.  According  to  Eq.  (4.1),  a  160  dB  signal  would  become 
triangular  in  about  9  wavelengths.  For  a  500  Hz  signal,  this  is  about  18  ft. 

When  the  wave  becomes  triangular,  the  energy  that  it  contains  is  reduced  by 
the  ratio  of  the  mean  square  of  a  triangle  to  a  sine  wave,  which  is  2/3  or  nearly 
2  dB.  In  addition,  the  500  Hz  component  of  the  triangular  wave  now  has  its 
amplitude  reduced  by  the  factor  2/rr,  representing  a  4  dB  loss  in  signal.  Thus, 
very  roughly,  we  can  say  that  nonlinearity  has  cost  us  4  dB  in  500  Hz  signal 
level.  Approximately  2  dB  of  this  is  lost  in  frequency  conversion  and  another  2 
dB  is  lost  in  dissipative  processes. 

Vibration  Simulation 

Vibration  testing  is  easily  the  most  common  form  of  environmental  testing 
for  dynamic  loads.  In  a  sense,  it  is  a  relatively  simple  form  of  testing— the  vibra¬ 
tion  output  of  a  shaker  is  generally  uniaxial  and  highly  controllable.  It  is  the 
nearest  thing  that  an  electrical  engineer  will  find  in  dynamics  to  the  one- 
dimensional  circuits  he  is  accustomed  to  studying. 
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direction  of  travel 
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Fig.  26.  Development  of  finite  amplitude  sine  wave  into  triangular 
wave. 


We  must  hasten  to  add,  however,  that  the  simplicity  of  vibration  simulation 
lies  entirely  in  the  test  equipment,  because  as  an  environment,  vibration  is  much 
more  complex  than  the  acoustic  environment  we  have  just  discussed. 

For  example,  a  structural  or  equipment  subassembly  like  the  truss-spacecraft 
system  in  Fig.  16  is  quite  intimately  attached  to  the  ring  frame-shroud  structure. 
The  dynamics  (or  impedances)  of  the  attachment  points  are  quite  important  in 
determining  energy  exchange.  The  dynamics  of  the  truss  system  are,  therefore, 
dependent  on  the  medium  (i.e.,  the  ring  frame)  to  which  it  is  attached.  This  is  in 
contrast  to  (he  dynamics  of  panels  in  energy  exchange  with  a  sound  field.  The 
air  loading  on  the  structure  does  not  greatly  alter  the  structural  dynamics. 

The  simulation  of  a  vibration  environment  as  currently  practiced  can  be  ex¬ 
emplified  by  continuing  to  take  the  truss-spacecraft  system  as  an  example.  Sup¬ 
pose  that  a  measurement  of  ring  frame  acceleration  in  a  direction  normal  to  the 
shroud  spectrum  was  obtained  from  a  similar  vehicle  that  had  the  spectrum 
shown  in  Fig.  10.  We  have  decided  that  curve  B  should  be  the  test  spectrum,  a 
shaker  is  attached  to  the  foot  of  one  of  the  trusses,  and  the  test  spectrum  is 
maintained  at  that  position.  If  the  motion  at  the  application  point  is  controlled, 
however,  the  natural  modes  of  the  truss  are  those  of  a  structure  with  the 
excitation  point  supported,  which  is  an  alteration  of  the  system.  We  cannot 
avoid  this  by  monitoring  the  force  in  the  service  vehicle  and  applying  a  given 
force  [88]  because  then  the  truss  modes  are  those  of  an  unsupported  or  free 
system. 

We  may  try  to  avoid  this  dilemma  by  attaching  the  trusses  to  a  fixture  (which 
may  be  similar  to  the  ring  frame)  and  shaking  the  fixture  with  a  controlled 
acceleration.  Unfortunately,  the  fixture  may  have  one  or  several  resonances  in 
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the  frequency  range  of  interest  that  change  the  acceleration  spectra.  These  may 
be  removed  electrically  by  employing  notch  filters  [89]  or  mechanically,  by 
stiffening  the  fixture.  In  either  case,  however,  the  dynamics  of  the  system  has 
been  substantially  altered  from  the  original  system,  with  attendant  complications 
in  mechanical  design  of  the  text  and  signal  processing. 

Scharton  has  pointed  out  recently  [90]  that  it  may  be  possible  to  avoid  some 
of  the  difficulties  of  the  test  fixture  by  replacing  it  with  a  light  and  flexible 
structure-a  multi-modal  test  fixture.  The  large  number  of  modes  in  the  fixture 
would  tend  to  smooth  out  the  acceleration  spectrum,  and  its  panel-like  character 
would  tend  to  simulate  the  structural  interactions  between  the  truss  and  ring 
frame-shroud  more  accurately.  This  approach  appears  to  have  considerable 
merit,  but  more  work  needs  to  be  done  to  prove  its  feasibility  in  general  test 
procedures. 

Aerodynamic  Simulation 

The  testing  of  a  vehicle  for  aerodynamic  loads  (in  distinction  to  determina¬ 
tion  of  loads)  is  difficult  to  do  by  direct  simulation.  If  simulation  is  attempted 
in  a  wind  tunnel,  several  problems  occur.  First,  very  few  wind  tunnels  are  large 
enough  to  accommodate  vehicles  of  interest.  Second,  the  large  amount  of 
acoustic  noise  existing  in  most  tunnels  interferes  with  the  interpretation  of 
response  to  aerodynamic  loads.  Finally,  wind  tunnel  operation  and  logistics  are 
very  expensive  and  complicated  compared  to  the  acoustic  and  vibration  facil¬ 
ities  we  have  discussed. 

In  at  least  one  case,  aerodynamic  simulation  of  loads  has  been  shown  to  be 
feasible.  In  a  study  of  the  excitation  of  a  flight  vehicle  by  a  turbulent  boundary 
layer,  Lyon  has  shown  that  it  should  be  possible  to  simulate  the  resonant  excita¬ 
tion  of  the  structure  by  a  suitable  distribution  of  normally  impinging  jets,  called 
wall  jets.  Each  wall  jet,  in  effect,  behaves  as  a  small  broadband  shaker.  These  are 
distributed  over  the  structure  in  such  a  way  that  a  fairly  homogeneous  vibration 
field  is  obtained.  The  spectrum  is  controlled  by  the  nozzle  diameter  and  volume 
flow  of  the  jets.  Details  are  given  in  Ref.  91. 

4.4  Simulation  by  Environment  Substitution 

It  often  happens  that  one  wishes  to  perform  an  environmental  test  using  an 
environment  different  from  the  one  being  simulated.  We  might  wish  to  avoid 
the  difficulties  just  cited  of  aerodynamic  simulation  by  exciting  the  structure 
with  an  equivalent  sound  field.  Or,  one  may  wish  to  replace  an  acoustic  excita¬ 
tion  by  a  source  of  mechanical  vibration,  cither  because  one  does  not  have  the 
required  acoustic  facility  or  because  the  vibration  test  is  judged  more  convenient. 
In  this  section,  we  explore  the  work  that  has  gone  on  in  substitution  of  loads,  in 
particular  examining  the  sense  in  which  one  load  is  equivalent  to  another. 
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Substitution  of  Acoustic  for  Aerodynamic  Loads 

The  aerodynamic  loads  that  one  wishes  to  simulate  are  usually  a  TBL, 
separated  flow,  or  oscillating  shock.  The  spatial  correlations  for  these  fields  are 
different  from  each  other,  and  they  are  also  different  from  that  of  a  sound  field. 
Thus,  in  general,  to  simulate  the  resonant  structural  response  due  to  the  TBL, 
we  would  require  a  sound  field  that  had  a  different  frenuency  spectrum  from 
that  of  the  TBL.  We  might  also  be  interested  in  the  sound  generated  within  the 
structure  by  the  TBL.  Since  the  internal  noise  depends  heavily  on  the  trans¬ 
mission  of  energy  by  nonresonant  modes  (recall  the  discussion  of  shroud  trans¬ 
mission  in  chapter  3),  the  differing  spatial  correlations  mean  that  if  the  frequency 
spectra  are  adjusted  to  give  the  same  resonant  response  spectrum  for  the  two 
loads,  then  the  nonresonant  tram  --on  will  likely  not  be  the  same. 

The  establishment  of  an  equivalent  sound  field  for  a  given  aerodynamic  load 
will  depend  on  the  kind  of  response  e  are  interested  in  and  the  nature  of  the 
structure.  The  equivalence  may  be  established  theoretically,  as  is  done  by 
Chandirami  et  al.  [2]  for  launch  vehicles,  or  it  may  be  done  empirically.  To  do 
the  latter,  one  requires  two  pieces  of  information:  (a)  field  data  which  relate 
TBL  loads  to  the  desired  response,  and  (b)  a  second  (perhaps  laboratory)  study 
which  relates  an  external  sound  field  to  that  same  response.  Neither  of  these 
fields  need  simulate  the  expected  environment  in  amplitude  or  spectral  shape, 
although  clearly  the  same  frequency  range  must  be  covered.  Once  the  expected 
TBL  spectrum  is  known,  it  is  converted  to  a  desired  response  spectrum.  This  ex¬ 
pected  response  spectrum  is  then  converted  to  an  acoustic  test  spectrum  using  the 
acoustic  excitation  response  transfer  function  that  was  experimentally  developed. 

As  an  example  of  the  foregoing,  suppose  we  are  asked  to  test  the  spacecraft 
shown  in  Fig.  16  for  its  response  during  the  transonic  portion  of  its  flight  profile. 
Figure  27  shows  the  acoustic  levels  measured  external  to  and  within  the  shroud 
at  launch  [92] .  The  difference  in  these  is  the  acoustic  noise  reduction  (NR)  of 
the  shroud  and  is  shown  in  Fig.  28. 

In  Fig.  29,  we  show  the  sound  levels  measured  within  the  shroud  curing  the 
shroud  during  the  transonic  portion  of  flight.  If  the  spacecraft  is  the  object  to 
be  tested,  then  these  levels  alone  would  be  sufficient  for  testing  purposes,  since 
they  could  be  applied  directly  to  the  payload.  If  the  overall  sluoud-spacecraft 
system  is  of  concern,  however,  then  the  NR  previously  obtained  should  be  added 
to  the  internal  levels  to  generate  an  equivalent  acoustic  test  spectrum  for  the 
transonic  flight  regime.  This  test  curve  is  also  shown  in  Fig.  29. 

Equivalence  between  Sound  and  Vibration 

There  are  many  reasons  why  one  might  wish  to  replace  a  specified  acoustic 
test  with  a  vibration  test,  and  vice  versa.  The  one  most  frequently  cited  is 
availability  of  test  facilities,  but  cost,  convenience,  and  ability  to  analyze  the 
results  are  also  considerations.  Many  aspects  of  the  substitution  problem  are 
still  in  the  research  phase.  Consider  cost,  for  example.  It  is  probably  less 
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Fig.  27.  Sound  pressure  levels  measured  external  to  and  within  a  spacecraft  shroud. 
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Fig.  28.  Acoustic  noise  reduction  derived  from  data  in  Fig.  27. 


expensive  to  generate  103  watts  of  power  into  an  acoustic  wave  than  it  is  to 
develop  the  same  mechanical  power  from  a  shaker.  In  an  environmental  tr 
however,  only  a  small  fraction  of  the  acoustic  power  is  absorbcd'by  the  struc¬ 
ture  from  the  sound  field.  Thus,  the  cost  per  watt  of  power  into  the  structure 
may  well  be  less  from  the  shaker. 
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The  most  comprehensive  effort  to  date  on  the  replacement,  of  acoustic  excita¬ 
tion  by  mechanical  shakers  is  the  work  of  Noiseux  [10] .  By  running  a  series  of 
laboratory  tests  on  the  airborne  computer  shown  in  Fig.  30,  Noiseux  was  able  to 
show  that  in  certain  frequency  bands  the  excitation  of  circuit  boards  was  due  to 
the  acoustic  environment,  and  in  others  it  was  due  to  the  vibration  at  the 
mounting  feet  positions  shown  in  Fig.  30.  The  combined  environment  was  then 
simulated  in  a  vibration  test  by  requiring  that  the  circuit  board  have  its  antic¬ 
ipated  vibration  levels.  A  comparison  between  realized  and  expected  vibration 
levels  on  the  circuit  board  is  shown  in  Fig.  31. 

A  major  difference  between  acoustic  and  mechanical  excitation  of  extended 
structures  is  the  greater  spatial  homogeneity  of  the  acoustic  excitation.  To  model 
acoustic  excitation  with  vibration,  therefore,  we  must  (a)  use  many  shakers, 
distributed  over  the  structure;  (b)  excite  an  intermediate  structure  (a  fixture) 
with  the  shaker  and  let  it  distribute  the  mechanical  energy  over  the  test  struc¬ 
ture;  (c)  rely  on  the  reverberation  to  diffuse  the  vibrational  energy  over  the 
structure,  or  (d)  some  combination  of  these. 

As  an  example  of  this,  suppose  we  are  interested  in  the  transmission  of  vibra¬ 
tion  to  the  spacecraft  shown  in  Fig.  16  along  its  mounting  truss  due  to  acoustic 
excitation  of  the  shroud.  We  may,  however,  elect  to  excite  the  shroud  with  a 
shaker  rather  than  a  sound  field,  in  part  to  minimize  the  acoustic  transmission  to 
spacecraft  through  the  interior  volume.  (A  shaker  will  produce  less  sound 
..cnismission  due  to  nonresonant  vibration  than  will  the  sound  field  [7,  Fig.  14]. 
We  may  want  to  give  ourselves  some  insurance  in  this  regard  also  by  placing  some 
absorption  in  the  interior  volume.) 
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Fig.  31.  Comparison  of  circuit  board  vibrations  predicted  from  environment  and  response 
studies  with  levels  measured  during  test  firing  [  10] . 

There  may  be  some  inhomogeneity  in  the  vibration  field  near  the  shaker,  due 
to  its  direct  field.  The  extent  and  magnitude  of  this  direct  field  is  dependent  on 
the  damping  of  the  shroud  and  on  the  internal  mechanical  impedance  of  the 
shaker  as  compared  to  the  point  input  impedance  of  the  shroud.  These  factors 
are  discussed  very  thoroughly  in  Ref.  94.  The  important  fact  for  the  present 
discussion  is  that  the  pomt  of  attachment  of  the  shaker(s)  to  the  shroud  should 
be  sufficiently  removed  from  the  ring  frame  so  that  ihe  ring  frame  is  exposed  to 
a  fairly  homogeneous  reverberant  vibrational  field  in  direct  simulation  of  its 
condition  under  acoustic  and  aerodynamic  excitation. 

Closure 

There  are  several  encouraging  features  of  the  present  patterns  of  development 
in  sound  and  vibration  environmental  studies.  One  is  the  attempt  to  gather  much 
more  environmental  data.  Such  data  are  essential,  both  for  developing  empirical 
and  semi-empirical  prediction  schemes  and  for  correlation  with  predictions  using 
rational  schemes  of  estimation. 

A  second  feature  that  is  encouraging  is  the  attempt  by  environmental  testing 
groups  to  develop  and  use  their  tests  for  evaluation  of  excitation  and  response 
parameters,  i.c.,  to  engage  in  relatively  more  analytic  testing  as  compared  to 
simple  proof  testing.  This  is  in  part  motivated  by  a  desire  to  reduce  or  modify 
test  specifications,  but  it  represents  a  desirable  trend  toward  seeking  better  under¬ 
standing  of  the  system  dynamics. 
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There  are  other  encouraging  features,  but  one  that  stands  out  particularly  is 
that  some  university  departments  are  now  training  students  in  the  fields  of 
aerodynamic  noise,  random  vibration,  and  the  response  of  complex  structures 
to  sound  and  mechanical  excitation.  These  fields  have  been  badly  neglected  in 
schools,  but  we  can  hope  that  many  more  will  develop  activities  in  environ¬ 
mental  sound  and  vibration  problems. 
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DEFINITION  OF  SYMBOLS 


A  Intercept  of  pressure-acceleration  level  correlation 

At  Turbulence  correlation  area 

a  Skin  acceleration 

a  Vehicle  radius 

BPF  Base  pressure  fluctuations 

Cprms  Fluctuation  pressure  coefficient 

c  Sound  speed  in  medium 

ct  Longitudinal  wave  speed  in  structural  material 

D  Jet  nozzle  diameter 

D  PSD  of  white  noise  process 

Ea  Total  energy  of  system  a 

f  Frequency 

f  Resonance  frequency  of  resonator 

f  Center  frequency  of  oscillatory  shock  motion 

f  Ring  frequency  of  cylindrical  vehicle 

hM  Modal  impulse  response 

Imn  Joint  acceptance  of  modes 

K  Scale  factor 

k  Acoustic  wave  number 

k  Boltzmann’s  constant 

ki  Wave  number,  conjugate  to  ,v1 

L  Distance  from  downstream  vehicle  flare  to  oscillating  shock 

L  Elastic  restoring  force  operator 

L0  Acceleration  level  (dB  re  1  g) 

Lp  Fluctuating  pressure  level  (dB  rc  0.0002  jubar) 

Lw  Acoustic  power  level  (dB  re  10"12  watt) 

Li  Correlation  scale  in  cross-flow  direction 

M  Mach  number  U/c 

M  Resonator  mass 

M  Slope  of  pressure-acceleration  level  correlation 

Ma  Mass  of  system  a 

n  Modal  density 

P(  )  Pressure  cross-PSD 
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Find  Pressure  in  induction  field 

PM  Pressure  amplitude  in  modal  expansion 

prej  Reference  pressure,  0.0002  pbar 

prms  Rms  of  fluctuating  pressure 

PSD  Power  spectral  density 

p(x,t )  Fluctuating  pressure 

q  Dynamic  pressure,  (l/2)pU2 

Rj y  Reynolds  number  per  foot 

Rrad  Radiation  resistance 

r  Radial  coordinate 

SEA  Statistical  energy  analysis 

S(f )  PSD  of  random  force  applied  to  resonator 

Sp(f)  PSD  of  pressure  excitation 

T  Absolute  temperature 

TBL  Turbulent  boundary  layer 

t  Temporal  coordinate 

U  Jet  efflux  speed 

Uc  Convection  speed 

V  Acoustic  volume 

V  Forwaid  speed  of  vehicle 

W  Structural  surface  density  (lb/sq  ft) 

WM  Displacement  amplitude  in  modal  expansion 

w  System  displacement 

x  Spatial  coordinate 

xc  Distance  to  shock  formed  by  high-intensity  sound  wave 

yM  Modal  admittance  function 

a  Angle  of  attack  of  vehicle 

0  Coupling  factor  relating  pressure  to  modal  force  in  Eldred’ 

procedure 

0  Viscous  damping  parameter 

A  Frequency  bandwidth 

A P  Pressure  jump  across  shock 

6*  Boundary  layer  displacement  thickness 

6(  )  Dirac  delta  function 

Kronecker  delta 
tj  Damping  loss  factor 

Damping  loss  factor  of  system  a 
qa  b  Coupling  loss  factor 

0  Eddy  lifetime  in  turbulence 
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^ a 

Ms 

ndiss 

ll/i 


n rad 

P 

P 

T 

<Pp 

^3 

“>M 

M 


Acoustic  wavelength 

Longitudinal  wavelength  in  structure  material 

Spatial  separation  in  downstream  and  transverse  directions 

Power  dissipated  in  system  a 

Power  flow  from  system  a  to  system  b 

Power  flow  from  system  a  to  system  b  by  nonresonant  motion  of 
intervening  system 
Radiated  sound  power 
Fluid  density 
Structural  mass  density 
Time  delay 

“Moving  axis”  temporal  correlation 
Space-time  pressure  correlation 
Downstream  space-time  correlation 
Cross-stream  spatial  correlation 
Modal  shape  function 
Undamped  modal  resonance  frequency 
Damped  modal  resonance  frequency 


<. . .  >  Averaging  operator 
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